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CHAPTER  1 


INTRODUCTION 

In  this  chapter,  I  sketch  a  short  history  of  work  leading  to 
this  thesis.  Then  I  give  a  summary  of  the  results  and  organization  of  the 
thesis . 

1.1.  Background 

Narendra  and  Valavani  (1978)  presented  a  stable  model 
reference  adaptive  controller.  Since  then,  it  has  been  recognized  that  seme 
of  the  assumptions  about  the  plant  or  about  our  knowledge  of  the  plant  are 
unrealistic.  It  has  become  popular  to  study  model  reference  adaptive 
systems  which  are  formed  by  applying  this  Narendra-Valavani  controller  to 
a  plant  which  meets  these  unrealistic  assumptions  and  then  perturbing  the 
plant  without  changing  the  controller. 

In  this  way,  researchers  have  been  able  to  study  more  realistic 
applications  of  adaptive  control  without  having  to  develop  entirely  new 
controller  structures.  In  this  thesis,  I  join  the  group  of  researchers  who 
have  perturbed  the  plant  by  adding  high  frequency  dynamics.  The  results  of 
a  few  of  the  researchers  who  have  preceded  me  in  studying  the  effects  of 
this  perturbation  are  summarized  in  the  following  paragraphs. 

Rohrs  et  al.  (1981,1982)  studied  the  effects  of  an  additive 
disturbance  and  unmodelled  high  frequency  dynamics  when  the  model  reference 
adaptive  system  is  excited  by  a  constant  reference  input.  They  show  that 
the  magnitude  of  the  adjustable  feedback  gains  in  the  controller  will 


eventually  become  large.  If  the  relative  degree  of  the  plant  including  the 


high  frequency  dynamics  is  greater  than  two,  the  linear  system  with  feedback, 
gains  fixed  at  some  large  values  would  be  unstable.  The  behavior  of  the 
model  reference  adaptive  systems  studied  by  these  researchers  is  characterized 
by  two  different  intervals  of  time.  During  the  first  interval  the  output 
error  becomes  and  remains  small  while  the  magnitudes  of  the  adjustable 
feedback  gains  drift  toward  infinity.  At  some  instant  the  adjustable  gains 
take  on  values  at  which  the  linear  system  with  feedback  gains  held  constant 
at  these  values  becomes  unstable.  This  is  the  beginning  of  the  second  interval 
during  which  the  magnitude  of  output  error  and  the  magnitudes  of  the  adjustable 
feedback  gains  approach  infinity  in  finite  time. 

The  same  group  also  studied  by  simulation  the  effects  of  high 
frequency  dynamics  when  the  reference  input  contains  high  frequencies  or 
the  initial  conditions  on  the  system  states  are  large.  They  show  that  both 
of  these  conditions  can  lead  to  instability.  However,  these  instabilities 
are  not  as  easily  described  as  the  drift  instability. 

The  drift  behavior  of  the  parameters  in  the  presence  of  an 
additive  disturbance  does  not  depend  on  the  high  frequency  dynamics.  This 
drift  occurs  because  the  control  scheme  adjusts  the  feedback  gains  until 
the  output  error  is  zero.  For  a  constant  reference  input,  there  exists 
an  unbounded  manifold  on  which  the  output  error  is  zero  when  the  disturbance 
is  zero.  Except  for  certain  special  disturbances,  perfect  disturbance 
rejection  requires  infinite  feedback  gain.  Hence,  the  feedback  gains  drift 
to  infinity  along  the  manifold.  This  behavior  is  discussed  by  Egardt  (1980) 
and  also  by  Riedle,  Cyr,  and  Kokotovic  (1983). 
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One  obvious  way  to  avoid  the  drift  to  infinity  is  to  put  a 
dead  zone  in  the  adaptation  law  around  output  error  equal  to  zero.  One  such 
algorithm  was  proposed  by  Peterson  and  Narendra  (1982). 

Ioannou  and  Kokotovic  (1983)  propose  a  different  solution. 

They  replace  the  integrator  in  the  original  adaptation  laws  with  a  first 
order  low  pass  filter.  For  a  bounded  additive  disturbance  the  output  error 
will  enter  and  remain  in  a  small  segment  containing  zero  and  the  adjustable 
feedback  gains  will  remain  bounded.  For  a  perturbation  involving  high 
frequency  dynamics,  they  define  the  term  "dominantly  rich"  input.  This 
is  sufficiently  rich  input  which  does  not  excite  the  high  frequency  dynamics. 
Then  they  show  that,  when  the  reference  input  is  dominantly  rich  and  the  system 
is  initially  in  a  region  of  attraction,  the  output  error  enters  a  small  segment 
of  the  real  line  containing  zero  and  the  adjustable  gains  remain  bounded. 

Krause  (1983)  refines  the  notion  of  dominant  richness  and 
studies  the  effect  of  high  frequency  dynamics  when  the  reference  input  is 
periodic,  has  at  least  as  many  spectral  components  as  unknowns,  and  is 
dominantly  rich.  He  finds  that  when  the  states  of  the  system  and  the 
errors  in  the  values  of  the  adjustable  gains  are  0(1) ,  the  adaptation  will 
be  in  the  correct  direction  if  the  speed  of  adaptation  is  slow  enough. 

The  approach  taken  in  this  thesis  to  make  the  system  robust  with 
respect  to  perturbations  is  to  make  the  unperturbed  system  exponentially  stable 
Kosut  (1983);  Kosut,  Johnson,  and  Anderson  (1983);  and  Anderson  and  Johnstone 
(1981)  have  also  used  this  approach.  Kosut  et  al.  introduced  the  concepts  of 
"tuned  system"  and  "tuned  error"  which  are  important  in  this  thesis.  They 
showed  that  persistent  excitation  which  provides  exponential  stability  of  the 
unperturbed  system  provides  robustness  with  respect  to  perturbations  with 
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bounded  effects.  However,  they  did  not  specify  any  particular  type  of 
perturbation  and  their  method  of  proof  is  considerably  different  from 
that  of  this  thesis.  Anderson  and  Johnstone  considered  the  discrete  time 
case  and  used  converse  Lyapunov  function  results  for  difference  equations  in 
a  manner  similar  to  the  way  I  use  converse  Lyapunov  function  results  for 
differential  equations. 

While  none  of  the  results  of  this  thesis  are  new,  I  feel  that 
this  work  represents  a  good  combination  of  the  techniques  and  ideas  mentioned 
above.  Choosing  a  particular  perturbation,  high  frequency  unmodelled 
dynamics,  allows  me  to  study  the  effects  of  the  perturbation  in  some  detail. 
With  the  concepts  of  tuned  system  and  tuned  error  I  am  able  to  specify  and 
solve  a  homogeneous  problem  before  attacking  the  complete  problem.  Using 
the  singular  perturbation  approach  of  Ioannou  and  Kokotovic  to  represent  the 
unmodelled  high  frequency  dynamics  provides  a  natural  way  to  partition  the 
system.  This  natural  partitioning  of  the  system  can  be  combined  with  the 
converse  Lyapunov  results  to  get  Lyapunov  functions  for  each  subsystem.  The 
final  idea  is  then  to  use  a  composite  of  the  Lyapunov  functions  for  the 
subsystems  to  prove  the  desired  stability  results. 

1.2.  Summary  of  Results  and  Organization  of  Thesis 

This  thesis  contains  two  major  results,  as  well  as  illustrative 
simulation  examples.  The  first  major  result  is  that  for  a  certain  class  of 
inputs  there  exists  a  unique  equilibrium  for  a  model  reference  adaptive 
system  when  the  nominal  plant  is  perturbed  by  high  frequency  unmodelled 
dynamics.  The  second  major  result  is  that,  for  a  small  enough  perturbation, 
this  equilibrium  is  exponentially  stable. 


In  Chapter  2,  I  derive  the  first  major  result  in  Theorem  2.1. 
This  chapter  also* provides  a  transfer  function  description  of  the  system 
and  contains  a  section  showing  that  a  small  singular  perturbation  has  a 
small  effect  on  the  equilibrium  values  of  the  adjustable  gains.  Chapter  2 
ends  with  remarks  which  are  intended  to  provide  either  an  intuitive 
interpretation  of  the  results  on  the  chapter  or  a  comparison  with  the  work 
of  others. 

In  Chapter  3,  I  derive  several  results  including  the  second 
major  result  of  the  thesis.  The  chapter  begins  with  definitions  and 
theorems  to  be  used  in  the  derivations  in  the  rest  of  the  chapter.  Then 
the  differential  equations  describing  the  system  are  presented  and  an 
error  system  is  derived.  When  the  tuned  error  is  zero,  as  it  will  be  for 
the  equilibrium  of  Theorem  2.1,  the  second  major  result  of  the  thesis  is 
stated  in  Theorem  3.9.  After  the  nonzero  tuned  error  case  is  handled. 
Chapter  3  ends  with  remarks  interpreting  the  results  of  the  chapter. 

Cyt  '  ; 

In  Chapter  A,  I  present  two  simple  exaraplesr which  show  that 
the  system  indeed  remains  exponentially  stable  when  high  frequency 
unmodelled  dynamics  perturb  the  original  system.  The  examples  also  show 
that  estimates  of  the  range  of  stable  perturbations  based  upon  the  proof  of 
Theorem  3.9  are  so  conservative  that  these  estimates  are  not  of  practical 
use . 

In  Chapter  5,  I  offer  some  concluding  remarks  and  suggestions 


for  future  research. 


CHAPTER  2 


EXISTENCE  OF  AN  EQUILIBRIUM 

In  this  chapter  I  have  three  major  objectives.  First,  I 
describe  the  perturbation  from  the  ideal  which  I  have  studied.  Next,  I 
show  that  for  a  certain  class  of  input  signals  a  unique  equilibrium  point 
exists  in  spite  of  the  perturbation.  Finally,  I  discuss  the  effects  of 
the  perturbation  on  the  parameter  values  associated  with  the  equilibrium. 

2.1.  Preliminaries 

The  following  result  will  be  useful  in  the  proof  of  the 
uniqueness  of  the  equilibrium  point. 

Lemma  1 :  Let  T^(s)  be  a  polynomial  of  degree  <  n-2 ,  T9(s)  be  a  polynomial  of 
degree  <  n,  R(s)  be  a  polynomial  of  degree  n,  and  Z(s)  be  a  polynomial  of  degree 
n-1.  Furthermore,  let  Z(s)  and  R(s)  be  relatively  prime.  Define  Q(s)  as  a 
polynomial  of  degree  <  2n-l,  such  that 

Q (s)  =  T1(s)R(s)  +  T2(s)Z(s)  .  (2.1) 

Let  X(us)  be  a  function  which  can  be  represented  by 

°o 

X(us)=l+  Z  x.(us)1Vs£{s:|s|<— }  .  (2.2) 

i=l  1  u 

Further,  require  that  X(ps)  satisfy 

-  x.  u  !  s|  <  I  x(ps)  I  <1  +  x„y I s  j  Vs  £  (s  :  I s I  <  — } 

li'i  1  —  )  1  1  1  1  II 


0  <  1 


(2.3) 


Define 


P(u,s)  =  T1(s)R(s)  +  T?(s)Z(s)X(ps)  .  (2.4) 

Then,  there  exists  UQ  Gl+,  such  that  ¥u€  [0,u  ),  if  P(u,s)  =  0  for  2n 
distinct  values  of  s;  s^Sc  i=l ,2 , . . .  ,2n,  then  T^(s)  =  0  and  T2(s)  =  0. 

Proof ;  First  consider  u  =  0.  Then  P(0,s)  =  Q(s).  Since  the  degree  of  Q(s) 

<  2n-l  and  Q(s)  is  zero  for  2n  distinct  values  of  s,  then  Q(s)  must  be 
identically  zero.  Therefore, 


T1(s)R(s)  =  -  T2(s)Z(s)  Vs  6  C 


(2.5) 


Since  R(s)  and  Z(s)  are  relatively  prime  and  neither  is  identically  zero, 
(2.5)  implies  either  T^s)  =  0  and  T2(s)  =  0,  or  T-^s)  contains  Z(s)  as  a 
factor.  The  latter  is  clearly  impossible  because  of  the  degree  of  Z(s)  is 
n-1  and  the  degree  of  T^(s)  is  strictly  less  than  n-1.  Thus,  the  lemma  is 
proved  for  u  =  0. 

Now,  consider  u>0.  Assume  T^(s)  i  0  or  T9(s)  t  0.  Using 
(2.1),  (2.3)  and  (2.4)  write 


P(u,s)  =  Q(s)  + 


T2(s)Z(s) 


(\x. 

W 1 


(us)  |  Vs  e  { s :  |  s 


(2.6) 


By  assumption  Q(s)  has  degree  greater  than  n-2  and  less  than  2n,  and,  hence, 
there  is  at  least  one  s^  such  that  Q(s^)^Q.  Without  loss  of  generality  let 
this  be  s^.  Let 


uo  =  min  {*T7 


Iq(s1) 


xl |T2(s1)Z(s1) 


’  T 


■  r  . 

a" 


si 


} 


(2.7 


Then,  if  u  6  (0,mq) 

|P(M,s1)|  >  |Q(s1) |  -  yx | | | T2 (s1) Z(s1) |  >  0  ;  (2.8) 

this  contradicts  P(y,s^)=C.  Therefore,  the  assumption  that  T1(s)j(0  or  T^CsJi^O 
is  false  and  the  lemma  is  proved. 


2.2.  The  System 

The  Model  Reference  Adaptive  System  to  be  studied  is  the 
Narendra,  Valavani  (1978)  controller  applied  to  a  plant  with  unmodelled  high- 
frequency  dynamics.  This  system  has  also  been  studied  by  Ioannou  and 
Kokotovic  (1983),  Rohrs  et  al.  (1981,1982),  and  Krause  (1983).  The  system 
has  the  block  diagram  shown  in  Figure  2.1. 

The  linear  time  invariant  plant  can  be  represented  by  a  series 
connection  of  two  transfer  functions,  Wp(s)  and  W^(us).  The  transfer  function 
of  the  nominal  or  slow  part  of  the  plant  can  be  written 


Z_(s) 

W  (s)  =  k  -~r — r* 
P  p  R  (s) 


(2.9) 


where  k^  is  a  scalar,  and  Z^(s)  and  Rp(s)  are  relatively  prime  monic  polynomials 
of  degree  n-1  and  n,  respectively.  The  transfer  function  of  the  fast  dynamics 


(2.12) 


Rf (ms)  =  1  +  Z  rf,  (ms) 
k-1 


The  following  assumptions  are  made  about  the  plant: 
Al)  n  is  known, 


A2)  the  sign  of  is  known  (without  loss  of  generality  let 
it  be  +) , 

A3)  Zp(s)  is  Hurwitz, 

A4)  |  Zf  (ms)  |  >  1  -  z1M  |  >0  Vs  £  {  s  :  |  s  |  <  -jj}  , 

A5)  if  Rc(us)  =  0  then  Re(s)  < -  —  , 

A6)  N2  -  1  ’  rf(Sl  +  s2)  *  °  ■ 

The  model  transfer  function  can  be  written 


Z  (s) 

VS>  "  km  FW 

m 


(2.13) 


where  k  is  a  positive  scalar,  Z  (s)  and  R  (s)  are  monic  Hurwitz  polynmials 
m  mm 

of  degree  n-1  and  n,  respectively.  Z  (s)  and  R  (s)  are  chosen  so  that  W  (s) 

'mm  m 

is  a  strictly  positive  real  transfer  function.  The  controller  has  the  structure 

T  T  T  9n 

shown  in  Figure  2.1  with  the  adjustable  parameter  vector  [cq,c  >dQ,d  J  e  1R~  . 

Z  (s)  in  the  controller  is  the  numerator  of  the  model  transfer  function  and 
m 


T(s)  is  defined 


T(s)  = 


.  n-2 
s 


(2.14) 


i 


'  ■  iifrfl  i 
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2.3.  The  Equilibrium 

In  this  section,  I  show  that, when  the  perturbation  parameter 
is  sufficiently  small  and  the  reference  input  r  belongs  to  a  certain  class 
of  signals,  an  equilibrium  exists  for  the  model  reference  adaptive  system 
described  in  the  previous  section.  As  usual,  let  the  output  error  be  the 
difference  between  the  plant  output  and  the  model  output,  e  =  y  -  y  . 

If,  for  some  constant  value  of  the  parameter  vector,  the  steady  state  output 
error  is  zero,  then  an  equilibrium  exists.  Hence,  when  an  equilibrium  exists, 
the  Fourier  transform  of  the  output  error  will  be  zero. 

Let  W(p,s)  be  the  transfer  function  of  the  controlled  plant 
from  the  reference  input  r  to  the  output  y  for  a  constant  value  of  the 
parameter  vector, 


W(y,s)  =  c  k 


Zm(s)Zp(s)Wf (us) 


°  P  Zm(s)Rp(s)  -cTt(s)rp(s>  -kp(doZm(s)+dTT(s))Zp(s)Wf(ps) 


.  (2.15) 


If  the  only  restrictions  on  the  reference  input,  r(t),  are  boundedness  and 
continuity,  then  for  the  steady  state  output  error  to  be  zero,  it  is 
necessary  that 


1 


W(p,s)  =  W  (s)  Vs  €  c  .  (2.16) 

m 

However,  when  P^O  the  relative  degree  of  W(p,s)  is  not  equal  to  the  relative 
degree  of  Wm(s) .  Therefore,  no  constant  parameter  vector  could  assure  (2.16). 

Clearly,  the  reference  input  will  have  to  belong  to  a  more 
restrictive  class  than  bounded,  continuous  functions  of  time.  Consider  the 
class  of  signals  which  are  sums  of  n  distinct  sinusoids,  that  is, 


(2.17) 


r(t)  =  l  ri  sin(wit  +  6i)  ,  r^O.w^O.w^w^d'-O,  ri,wi,0i 


e  k 


Denote  the  set  of  s-values  in  (2.17)  by  S,  that  is, 


S  =  {jw1,-jw1, . . . ,jwn,-jwn) 


(2.18) 


For  an  equilibrium  to  exist  with  an  input  signal  in  this  class,  it  is 
necessary  and  sufficient  that 


W(y,s)  =  W  (s)  for  sSs 
m 


(2.19) 


This  is  equivalent  to 


Zm(s)Rp(s)  =  c  T(s)Rp(s) +  kp 


t—  R  (s) +  d  Z  (s) +  d  T(s) 
k  m  o  m 


l  m 


Z  (s)Wf(ys), 


(2.20) 


s  €  S 


Theorem  2.1:  There  exists  y  €  ]R  such  that  ¥y  £  [0,Uq),  the  MRAS  shown 
in  Figure  1  and  described  by  (2. 9) -(2. 14)  with  assumptions  A1-A6  has  a 
unique  equilibrium  for  a  reference  input  described  by  (2.17). 

Proof :  (2.20)  is  a  set  of  2n  linear  equations  in  the  2n  elements  of  the 

parameter  vector.  Clearly,  there  exists  at  least  one  constant  parameter  vector 
which  satisfies  (2.20).  Furthermore,  if  the  2n  equations  (2.20)  are  linearly 
independent,  the  constant  parameter  vector  which  satisfies  (2.20)  is  unique. 

Assume  that  the  2n  equations  are  not  linearly  independent. 

Then  there  must  exist  a  nonzero  constant  parameter  vector  such  that 


0  =  c  T(s)R  (s) +  k 


-T-  R  (s)  +  d  Z  (s)  +dTT(s) 


Z  (s)W  (ps)  ,  se  S  .  (2.21) 


By  Lemma  1  there  exists  a  uq6]R  such  that,  Vy  G  [0,Lo),  the  only  solution 
to  (2.21)  is  the  zero  vector.  Thus,  Vy G  [0,yQ),  the  2n  equations  must  be 
linearly  independent  and  the  theorem  is  proved. 


2.4.  Effects  of  the  Perturbation 


Let  0  (0)  be  the  unique  constant  parameter  vector  for  which 


the  steady  state  error  is  zero  in  the  unperturbed  case,  y=0.  For 


T  T  T  * 

[cq,c  >dQ,d  ]  =  0  (0)  and  y=0,  the  transfer  function  of  the  controlled 


plant  is  identically  equal  to  the  transfer  function  of  the  model.  Since 


the  model  transfer  function  is  strictly  positive  real,  so  is  the  controlled 

k 


plant  transfer  function.  In  this  section,  I  want  to  show  that  8  (y) ,  the 
unique  constant  parameter  vector  for  which  the  steady  state  error  is  zero 
in  the  perturbed  plant,  satisfies 


0  (y)  =  6"(0)  +  0(y) 


(2.2 


k  k  'p  k 

I  also  want  to  show  that  for  [cq,c  ,dQ,d  ]  =  0  (u),  the  transfer  function 


* 

c  k 
o  p 


Z  (s)Z  (s) 
m  p _ 


Z  (s)R  (s)  -  c*TT(s)R  (s)  -  k  (d*Zm(s) +  d*TT(s) ) Z(s) 
m  p  p  p  o  m  p 


(2.2 


will  be  strictly  positive  real  provided  y  is  small  enough.  This  transfer 


function  is  the  transfer  function  of  the  controlled  unperturbed  n-th  order 


plant  when  the  parameter  vector  is  held  constant  at  0  (y) 


From  (2.20)  0  (p)  sarisfies 


zm(s)Rp(S> 


*  T 

6  (p) 


k  Rm(s)z  (s)wt(i.s) 
m 


T(s)Rp(s) 


kpzm(s)Zp(s)Wf(ys) 


kpT(s) Zp (s)Wf (us) 


s  S  S 


Define  the  2n  by  2n  matrix 


H(u)  = 


W  -  Rm(^wn)Zp(-jwn)Wf(-Mjwn) 

m  m 


l 


m 


T  ( j  wx  )  Rp  (j  wx) . T  ( j  wn)  Rp  ( -j  wn) 

Vm(jwl)Zp(jwl)Wf(lJ:iwl)*--*kPZni(":iwn)Zp('jwn)Wf(-  jwn> 

kpT(jWl)Zp(jWl)Wf (tJjWl) . kpT(-jwn)Zp(’jwn)Wf (~  Jwn) 


From  Theorem  2.1,  if  p£  [0,p  ),  H(p)  is  nonsingular.  Assume 


u  <  min  — } 

s  SS  II 


Then,  H(p)  has  a  series  representation, 


H  (p)  =  H(0)  +  E  H.u1 
i=l 


It  follows  from  (2.27)  that 

CO 

H(p) _1  =  H(0)_1  +  E  L.u1 


r 


Then, 


a  -IT 

e  (u)  =  (H(u)  ) 


Zm(jwl)Rp(jwl) 


Z  (-jw  )R  (-jw  ) 
m  n  p  n 


zm(jWi)Rp(jWi) 

00  .  T  I 

[H(0)_1  +  Z  L  u1]  =  : 

i=l  Z  (-jw  )R  (-jw  ) 

m  J  n  p  J  n 


0*(o)  +  z  e.u1 


(2.29) 


Thus,  the  desired  property  (2.22)  holds  for  p  sufficiently  small. 

For  a  transfer  function  h(s)  to  be  strictly  positive  real,  it 


must  sati.-fy: 


1)  h(s)  is  real  for  real  s, 

2)  the  poles  of  h(s)  should  lie  in  Re[s]  <  0, 

3)  for  all  real  w,  one  had  Re[h(jw)]  >0,  -  °°  <  w  <  °°. 


Property  (3)  is  equivalent  to 


3')  for  all  real  w,  one  has  -90°  <  Phase  [h(jw)]<90°  -  °°  <w<  <*>. 
Theorem  2.2:  There  exists  a  u  > 0  such  that,  if  p€[0,p^)  and 


*  *t  *  at  1  A 

[co,c  i,do,d  j  =  9  (p) 


(2.30) 


then  the  transfer  function  (2.23)  is  strictly  positive  real. 


Proof :  Condition  (1)  is  met  by  construction  of  (2.23).  The  poles  of  (2.23) 

are  the  zeros  of  the  polynomial 

Z  (s)R  (s)  -  c*TT(s)R  (s)  -  k  (d*Z  (s)  +  d*XT(s) ) (s) 
m  p  p  p  o  m  p 

0 

T(s)R  (s) 
k  Z  (s)Z  (s 

pm  p 

k  T(s)Z  (s) 

P  P 

The  zeros  of  (2.31)  are  continuous  functions  of  y .  Since  the  zeros  of 

R  (s)Z  (s)  lie  in  Re[s]<0,  there  exists  a  y'  such  that,  if  y€  [0,y'), 
m  p 

then  the  poles  of  (2.23)  lie  in  Re[s]  <  0.  The  zeros  of  (2.23)  are  not 
functions  of  y.  Since  the  phase  is  a  continuous  function  of  the  pole  positions, 
and  the  pole  positions  are  continuous  functions  of  y,  the  phase  is  also  a 
continuous  function  of  y.  The  transfer  function  (2.23)  satisfies  condition 
3'  for  y=o.  Hence,  there  exists  y"  >  0  such  that  for  y£  [0,y"  ),  condition 
3'  is  satisfied  for  the  transfer  function  (2.23).  Take  y^  =  min [y ' , y "  ]  and 
the  theorem  is  proved. 

2.5.  Discussion 

In  this  chapter,  I  have  shown  that  for  a  certain  class  of  inputs 
and  a  sufficiently  small  perturbation  in  the  structure  of  the  plant,  a  unique 
equilibrium  point  exists.  Furthermore,  it  is  possible  for  the  controlled 
nominal  plant  with  the  adjustable  parameters  held  constant  at  the  equilibrium 
values  to  be  strictly  positive  real.  This  second  fact  will  be  important  in 
determining  the  stability  of  this  equilibrium  point. 


(2.31) 


I  would  like  to  highlight  a  few  points  which  were  not  explicitly 


! 

I 


stated  in  the  derivation  of  this  result. 

*  * 

Remark  2.1:  The  perturbation  of  0  from  9  (0)  will  be  a  function  of  the 
frequencies  w^  in  S.  This  can  be  seen  from  (2.25).  It  can  be  intuitively 
seen  using  the  following  argument.  The  equilibrium  exists  because  the  Bode 
plot  of  the  controlled  plant  transfer  function  can  match  that  of  the  model  at 
the  n  frequencies  in  the  input  signal.  The  contribution  of  the  unmodelled 
dynamics  to  the  controlled  plant  Bode*  plot  is  a  function  of  frequency. 

Therefore,  the  amount  by  which  the  Bode  plot  of  the  controlled  nominal  plant 
is  moved  from  the  Bode  plot  of  the  model  is  a  function  of  the  frequency. 

It  is  not  hard  to  believe,  then,  that  choosing  different  values  of  w^  will 
result  in  different  perturbations  of  the  parameters. 

Remark  2.2:  Choosing  a  reference  input  with  more  than  n  distinct  frequencies 
will  result  in  a  system  that  has  no  equilibrium.  This  is  true  because  it 
can  be  shown  in  a  manner  very  similar  to  the  proof  of  Theorem  2.1  that  the 
resulting  set  of  more  than  2n  equations  (2.20)  in  2n  unknowns  is  inconsistent. 
One  can  also  think  of  this  in  terms  of  the  Bode  plots.  The  2n  parameters 
provide  2n  degrees  of  freedom  with  which  to  bend  the  Bode  plot  of  the  controlled 
plant.  It  is  not  possible  to  make  the  (2n  +  l)-th  bend  to  match  at  an  extra 
point . 

Remark  2.3:  While  the  theorems  in  this  chapter  prove  existence  of  uq,  pj  , 
they  can  be  used  in  a  different  manner.  For  a  given  p  which  is  small  enough 
so  that  the  nominal  plant  dynamics  and  the  unmodelled  dynamics  have  a 
reasonable  separation,  the  w.'s  of  the  reference  input  can  be  chosen  so  that 


the  unique  equilibrium  point  exists.  This  is  related  to  the  idea  of 
dominant  richness  which  was  introduced  by  Ioannou  and  Kokotovic  (1983). 
Krause  (1983)  discusses  how  to  assure  that  the  reference  input  is  in  the 
dominantly  rich  range. 


CHAPTER  3 


STABILITY  OF  THE  EQUILIBRIUM 

In  this  chapter,  I  study  the  stability  of  the  equilibrium  of 
Theorem  2.1  when  the  parameter  vector  is  allowed  to  be  time  varying.  After 
some  definitions  and  theorems  about  stability,  I  present  the  differential 
equations  which  describe  the  behavior  of  the  model  reference  adaptive  control 
system.  From  this  system  of  equations  I  derive  an  error  system  which  describes 
the  variations  about  the  equilibrium. 

The  error  system  will  have  an  input  which  is  the  difference  of 
the  model  output  and  the  controlled  plant  output  when  9=6*.  In  Chapter  2  it 
was  shown  that  for  an  input  with  exactly  n  sinusoids,  the  steady  state  value 
of  this  difference  will  be  zero.  In  this  chapter,  I  will  show  that  when 
this  difference  is  zero,  there  exists  a  u  such  that  for  all  fO,u),  the 
zero  solution  of  the  error  system  is  exponentially  stable.  Then  I  will  show 
that  this  difference  is  always  0(u)  and  that  trajectories  of  the  error  system 
which  begin  in  some  region  of  attraction  converge  exponentially  to  an 
0(a)  residual  set  containing  zero. 

3.1.  Preliminaries 

To  insure  that  the  reader  and  I  agree  on  terminology  of 
stability  concepts,  I  state  the  following  definitions.  These  definitions  are 
taken  from  Yoshizawa  (1966)  and  Rouche,  Habets,  Laloy  (1977).  Consider  a 
system  of  differential  equations 


Suppose  that  f(t,x)  is  continuous  on  IRx  ]R  .  Assume  that  (3.1)  has  solutions 
which  are  uniquely  determined  by  the  initial  condition  (t,x)=  (to,XQ).  Let 
x(t;xQ,to)  be  the  solution  of  (3.1)  which  passes  through  the  point  (to>XQ) . 
Definition  3.1:  The  solution  x(t)  = 0  of  (3.1)  is  uniform-stable,  if  for 
any  e  >  0  and  tQ  £  ]R,  there  exists  a  6(e)  >  0  such  chat  if  ||  x  ||  <  6(e),  then 

<  £  for  a11  t-to< 

Definition  3.2:  The  zero  solution  of  (3.1)  is  uniform-asymptotically  stable 
if  it  is  uniform-stable,  and  if  given  any  >  0  and  any  t  £  1R,  there  exist 
a  6^  >  0  and  a  T(e^)>0  such  that  if  |xq||  <6^,  then  ||x(t;x  ,tQ)|  <  e^  for  all 
tito  +  Kep. 

Definition  3.3:  The  zero  solution  of  (3.1)  is  exponentially  stable  in  the 
large,  if  there  exists  an  a  >  0  and  for  any  6  >  0  there  exists  a  K(3)  L  1 
such  that  if  ||  xj|  <  2 


-"('•vVII  sK(8)e"“<t"t°>KI 


(3.2) 


Definition  3.4:  A  function  a  :  R+-*  R+  is  said  to  be  an  element  of  class  K 
if  it  is  continuous,  strictly  increasing,  and  a(0)  =  0.  We  write  a£K. 

Let  V(t,x)  be  a  continuous  scalar  function  defined  on  lx  lRn 
and  let  V(t,x)  satisfy  locally  a  Lipschitz  condition  with  respect  to  x. 
Then  define  the  function 


D+,.  ,.V(t,x)  =  lim.  p-  { V  ( t  +  h ,  x  +  hf  (t  ,x) )  -  V(  t  ,x) } 


(3.3) 
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If  V(t,x)  has  continuous  partial  derivatives  of  the  first  order, 


D|3>l)v(t,x)  =  f^+ff  •  F(t,x)  , 


(3.4) 


where  denotes  scalar  product. 

The  following  theorems  will  be  used  to  investigate  the  stability 
of  the  equilibrium  of  Theorem  2.1.  The  theorems  are  stated  here  without  proof. 
Theorems  3.1,  3.2,  3.3  and  3.4  are  standard  results  and  proofs  can  be  found 
in  a  number  of  texts.  Theorem  3.5  is  a  well-known  result  of  Lasalle  and  the 
proof  can  be  found  in  Hale  (1969) .  Theorems  3.6  and  3.7  appear  in  a  monograph 
by  Yoshizawa  (1966).  The  proofs  of  Theorems  3.6  and  3.7  can  be  found  in  the 
appendix  to  this  thesis.  Consider  the  linear  system 


x  =  A(t)x 


(3.5) 


where  A(t)  is  an  nxn  matrix  of  continuous  functions  of  time  defined  on  1R. 
Theorem  3.1:  If  the  zero  solution  of  (3.5)  is  uniform-asymptotically  stable, 
then  it  is  exponentially  stable  in  the  large  and  K  in  (3.2)  can  be  chosen 
independent  of  3. 

Now  consider  the  differential  equation  (3.1)  under  the  assump¬ 
tions  that  f(t,x)  is  continuous  on  0  <  t  <  llxi!  <  H,  H>0  and  f(t,0)=0. 
Theorem  3.2:  Suppose  that  there  exists  a  Lyapunov  function  V(t,x)  defined 
on  0  <  t  <  “,  || x|j  <  H  which  satisfies  the  following  conditions; 

i)  V(t,o)  =  0, 

ii)  a  ( |  x|  )  <  V(t,x)  <  b(j|  xjj  )  ,  a€K  ,  bSK, 

iii)  d|3  1^V(t,x)  <  0. 

Then,  the  solution  x(t)  =  0  of  the  system  (3.1)  is  uniform-stable. 


B 
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Theorem  3.3:  Under  the  same  conditions  as  Theorem  3.2,  if 

D(3  ^VU.x)  <  -  c(|x||),  where  c^K 

then  the  solution  x(t)  3  0  of  (3.1)  is  unif orm-asymptotically  stable. 

Corollary  3.3:  Under  the  same  conditions  as  Theorem  3.2,  if 

D(3  1^)V(t,x)  <  -  cV(t,x)  , 

where  c>0  is  a  constant,  then  the  solution  x(t)  =  0  of  (3.1)  is  uniform- 
asymptotically  stable. 

Assuming  that  f(t,x)  is  continuous  on  0  <  t  <  <=°,  x€Rn,  and 
f(t,0)X0,  the  following  theorem  applies. 

Theorem  3.4:-  Suppose  there  exists  a  Lyapunov  function  V(t,x)  defined  for 
0  <  t  <  «,  x€  ]Rn  satisfying  the  following  conditions: 

i)  |  x  j|  <  V(t,x)  <  K(g)  |x|  for  |  xj|  <  6, 
ii)  ^V(t,x)  <  -  aV(t,x),  where  a>0  is  a  constant. 

Then,  the  zero  solution  of  (3.1)  is  exponentially  stable  in  the  large,  that  is, 
N^VVII  -  K(g)e"a(t"to)  JxJ  for|xJ<B. 

Theorem  3.5:  Let  V(t,x)  be  a  continuous  Lyapunov  function  defined  on 
0  <  t  <  00 ,  xSG,  where  G  is  an  open  set  in  Rn.  Suppose  that 

i)  given  x€  G,  the  closure  of  G,  there  is  a  neighborhood  of  x,  N 
such  that  V(t,x)  is  bounded  from  below  for  all  t>0  and  all  x  in  NxnG,  and 
ii)  d|3  1^V(t,x)  <  -  W(x)  <  0  for  0  <  t  <  ”,  x€G  and  W(x)  is 

continuous  on  G. 

Define 

E  =  {x:  x^G  and  W(x)  =  0}  .  (3.6) 


Let  x(t;x  ,t  )  be  a  solution  of  (3.1)  which  is  bounded  and  remains  in  G 
o  o 


for  t  >  t  >0. 

-  o  ~ 


If  W(x)  has  continuous  first  derivatives  on  G  and 


W •  f(t,x)  is  bounded  from  above  along  the  solution  x(t;x  ,t  ), 


x(t;x  ,t  )  -*■  E  as  t 
o  o 


Theorem  3.6:  Suppose  there  exists  a  K  >  1  and  a  c  such  that 


‘(“VVl  -  Ke'c<t'toV„l  • 


where  x(t;xo,tQ)  is  a  solution  of  (3.5)  and  c  is  a  constant  (>  0).  Then, 
there  exists  a  Lyapunov  function  V(t,x)  which  satisfies  the  following 
conditions : 

i)  j|x||  <  V(t,x)  <  K||x||  , 

ii)  |  V(t,x)  -  V(t,x")  |  <  K||x  -  x'j|, 

iii)  d|3  5^V(t,x)  <  -  cV(t,x). 

Theorem  3.7:  Suppose  that  f(t,x)  of  (3.1)  is  continuous  for  0  <  t  < 
x-Hnand  f(t,0)  =0.  If  f(t,x)  satisfies  a  local  Lipschitz  condition  with 
respect  to  x  and  the  zero  solution  of  (3.1)  is  exponentially  stable 
in  the  large,  i.e.,  there  exists  an  t  •  0  and  for  any  3  >  0,  there 


exists  a  K(3)  >  0  such  that  if  |x  |j  <8, 


then  there  exists  a  Lyapunov  function  V(t,x)  defined  for  0  <  t  <  <*>,  x€  ]Rn 
which  satisfies  the  following  conditions: 

i)  JxJ  <  V(t,x)  <  K(8)  |x|  for  JxJ  <  6, 

ii)  |  V(t  ,x)  -  V(t,x")  |  <  L(t,6)  I  x  -  x'J  for  |  x||  <  8,  ||x'||  <  3, 

iii)  .^V(t,x)  <  -  qaV(t,x),  where  0<q<l. 

3.2.  The  Error  System 

For  the  unperturbed  case,  M=0,  the  proof  of  stability  of  the 
& 

equilibrium  (e,0)=(O,0  )  uses  an  error  system  representation  of  the  model 
reference  adaptive  control  system.  In  this  section,  I  derive  an  error 
system  representation  for  the  perturbed  case.  It  is  a  specific  example  of  the 
type  of  error  system  representation  which  has  been  presented  recently  by 
Kosut ,  Johnson,  and  Anderson  (1983). 

My  derivation  will  be  made  in  several  steps.  First,  I  present 
the  differential  equations  describing  the  plant,  the  control  system,  and  the 
model.  Next,  I  use  a  singular  perturbation  type  of  transformation  on  the 
plant.  This  explicitly  separates  the  high  frequency  unmodelled  dynamics 
from  the  nominal,  or  slow  part,  of  the  plant.  The  third  step  is  the 
introduction  of  a  "tuned  system."  Finally,  I  replace  the  original  reference 
model  with  the  tuned  system  in  order  to  get  the  desired  error  system 


representation . 


The  differential  equations  describing  the  plant  are 


x  =  Aux  +  b1h2z  , 


uz  =  A222  +  t>2U  , 


xe  ir 


nx+n2 


y  =  h^x  . 


The  vectors  and  matrices  in  (3.8)— (3. 10)  satisfy  the  following  equations 

hi(sl  “  An)-lbi =  VS)  ’ 

h2(usl  -  A22)  1b2  =  Wf(ys) 


The  equations  describing  the  auxiliary  signal  generators  are 


■ 

.1 

V 

=  Av1  +  bu  , 

lp:  Dn~l 
v  €  R 

.2 

V 

2 

=  Av  +  by  , 

2  r-  n-1 
v  e  R 

■ 

where  the  vector  b  and  matrix  A  satisfy 

(si 

-1  _  T(s) 

Z  (S) 
m 

• 

!■* 

- 

Define  Che  parameter 

vector,  9,  and  the 

signal 

vector »  w , 

9  - 

T  t  ^ 

(c0-c  .d0.d  ] 

9 

-*■ 

w  = 

]T  oT  t 

[ r , v  ,v,v  ] 

• 

Ik 

The  parameter  adjustment  law  is  then 

9  =  -  Tw(v  -  v  ) 
'  m 


and  the  input  to  the  plant  is 


u 


(3.19) 


* 

Thus,  I  have  a  set  of  differential  equations  which  describe  the  behavior  of 
the  MRAS  shown  in  Figure  2.1. 

The  choice  of  a  representation  for  the  "unmodelled"  part  of 
the  plant  is  the  single  most  important  factor  in  a  study  of  the  effects  of 
"unmodelled"  dynamics.  I  have  restricted  my  unmodelled  dynamics  to  the  high 
frequency  range.  Because  of  this  I  was  able  to  choose  a  system  of  differential 
equations,  (3.8)— (3.9) ,  in  singular  perturbation  form  to  represent  the  plant. 
Because  there  are  stability  results  for  the  unperturbed  case,  I  would  like 
to  have  my  plant  represented  in  a  form  which  appears  as  an  unperturbed  plant 
plus  a  perturbation.  Singular  perturbation  theory  suggests  the  following 
method  to  transform  (3.8)— (3.9)  to  the  desired  representation.  First,  set 
p=0  in  (3.9).  This  corresponds  to  making  the  fast  part  of  the  plant 


infinitely  fast.  Because  k^  is  assumed  stable  z  converges  infinitely 
fast  to 

z  =  -A.~b„u 


Then  z  is  considered  Che  slow"  part  of  z. 


In  order  to  separate  the  system 


into  its  fast  and  slow  parts,  I  define  a  new  variable 


n  =  z  -  z  =  z  +  A22b2u 


(3.23) 


to  be  the  fast  variable.  Now  I  can  rewrite  (3.8)— (3. 9)  as 


Anx  +  b]u  +  b1h2n 


A22n  +  yA22V 


(3.24) 

(3.25) 


The  MRAS  using  (3 . 24) -(3 . 25)  has  a  block  diagram  shown  in  Figure  3.1. 

The  third  step  in  this  derivation  is  the  definition  of  a 
* 

"tuned  system."  Define  r  (t)  to  be  an  input  as  in  (2.17).  Assume  that 

*  *T  *  *T  T 

p€  [0,uq)  so  that  a  unique  equilibrium  exists.  Take  [c^.c  ,do>d  1  to 

* 

be  the  solution  of  (2.20)  with  r(t)=r  (t).  Then,  the  tuned  system  is  the 

linear  time  invariant  system  formed  by  the  plant  plus  the  controller  with 
X  X  _  *  *T  *  *T 

[cq,c  ,dQ,d  ]  =  [cq,c  >d0,d  1’  An  imPort;ant  point  is  that,  whenever  p^O, 

the  tuned  system  is  not  positive  real.  In  fact,  the  order  and  relative 
degree  of  the  tuned  system  are  unknown.  The  tuned  system  is  represented 
by  the  equations 


X  = 

Anx 

+  b^i 

.1* 

* 

V  = 

aV 

+  bu 

.2* 

2* 

•k 

V  = 

#  * 

Av 

+  by 

* 

n 

A„„n 

+  PA, 

.  T  * 

>lh2r 


-U  •* 


(3.26a) 
(3.26b) 
(3.26c) 
( 3 . 26d) 


w 


*  *T  * 

u  =  6  w  , 

*  T 

w  =  [r,v  ,y  ,v“  ] 


(3.28) 


(3.29) 


In  order  to  substitute  the  tuned  system  for  the  reference  model, 
I  must  introduce  the  concept  of  tuned  error.  The  tuned  error  is  denoted  by  e 
and  is  defined  as  the  difference  between  the  output  of  the  tuned  system  and 
the  reference  model. 


*  A  * 

e  =  y  -  y 


(3.30) 


Notice  that  the  tuned  error  is  defined  for  any  reference  input  r(t).  Also 

ft  ft 

notice  that  when  r(t)=r  (t) ,  the  steady  state  value  of  e  will  be  zero. 
Figure  3.2  illustrates  the  tuned  system  and  the  tuned  error. 

Now,  define  the  error  system  variables  as 


=  x  -  x  , 


(3.31a) 


1  1* 
v  -  v 


(3.31b) 


2  2* 
v  -  v 


(3.31c) 


=  e  -  e 


( 3 . 31d) 


=  n  -  n 


(3 . 31e) 


The  differential  equations  describing  the  behavior  of  the  error  system  are 


x'  =  A, , x '  +  b . u '  +  b,hnn' 


(3.32a) 


1  ’  1  ’ 
v  =  Av  +  buf 


(3.32b) 


.  2  1  2  ’  , 

v  =  Av  +  by  , 


(3.32c) 


r  w(y '  +  e  ) 


(3 . 32d) 


where 


A9?n'  +  uA??b2u’ 


.  J  ,  * 

y  =  h^x  =  y  -  y 

*  1  *  2 ' 
w'  =  w  -  w  =  [0,v  ,y '  ,v  ]  , 

*T  T 

u '  =  9  w '  +  w  6  '  , 


a'  3  w  9'  +  w  9'  +  9  w'  +  9'  w' 

'  ^  T 

=  -  w^Tw(y'  +  e  )  +  w  9' 

T  T  1  '  T  2  '  T  T  T  AT  T 
+  f[[x’  ,v  ,v  ,9’  ]  +  doh1b1h2n' 


T  1-T  o»T  x  t 

+  i-]Jf  *WiV 

2  i 

L  -j 

+  d’h^b  9,Tw’ 
oil 


(3.33) 


(3.34) 

(3.35) 

(3.36) 


(3.37) 


*  T  ,  ,  T  *  A  t  a  T  A  9 

d0Anhi  +  (b  d  +  dob  c  +  bi  b:V)bi 

,T  *  T  *  T  a  a 

A  c  +  (b  c  +  hx  b3do)c 

T  *  T  A  T  A 

A  d  +  (b  c  +  h:  b,_do) 

T  A  T  *  A 
(b  c  +h,b,d  )w 


(3.38) 


Two  differences  exist  between  the  error  sys  T  derived  for 
the  perturbed  case  and  the  error  system  for  the  unperturbed  case.  The  first 
is  the  perturbation  itself.  Because  of  the  perturbation,  the  linear  time 
invariant  part  of  the  error  model  is  not  positive  real.  The  second  difference 

k 

is  the  tuned  error,  e  ,  which  shows  up  as  an  extra  external  input.  In  this 

section  I  study  the  stability  of  the  zero  solution  of  the  error  system  when 

* 

e  =  0.  I  will  consider  the  effects  of  a  nonzero  tuned  error  in  Section  3.4. 

The  results  of  this  section  can  be  summarized  as  follows. 

Assume  that  e  H  0  and  the  reference  input  has  at  least  n  distinct  sinusoids 
in  the  low  frequency  range.  Then  the  zero  solution  of  the  error  system 

(3.32) ,  (3.33)  is  exponentially  stable.  Furthermore,  the  magnitude  of  the 

_  1 

region  of  attraction  for  this  exponentiallv  stable  equilibrium  is  0(p  2) 

_  3 

for  the  slow  variables  and  0(u  2)  for  the  fast  variables. 

These  results  are  derived  in  several  steps.  The  first  step 
is  to  decompose  the  error  system  into  a  fast  subsystem  and  a  slow  subsystem. 
The  decomposition  is  performed  so  that  the  fast  subsystem  is  a  stable  linear 
time  invariant  one  by  assumption.  The  next  step  is  to  prove  that  the  slow 
subsystem  is  exponentially  stable.  Then,  I  apply  Theorem  3.7  to  the  slow 
subsystem  and  Theorem  3.6  to  the  fast  subsystem  to  get  a  conceptual  Lyapunov 
function  for  each  subsystem.  Finally,  I  use  a  combination  of  these  two 
Lyapunov  functions  to  prove  that  the  zero  solution  of  the  error  system  (3.32), 

(3.33)  is  exponentially  stable  in  the  region  of  attraction  of  the  zero 
solut ion . 

I  get  a  good  decomposition  for  this  problem  using  the  techniques 
of  singular  perturbation  ^heory.  This  decomposition  was  set  up  by  the 


transformation  (3.23)  which  introduced  n  to  serve  as  a  fast  variable.  To  get 


the  slow  subsystem  from  the  error  system  (3.32),  (3.33),  I  set  u=0  in  (3.33) 
which  corresponds  to  making  my  unmodelled  dynamics  infinitely  fast.  Since 
A0O  is  stable  any  initial  conditions  on  n '  will  decay  to  zero  infinitely  fast. 
Hence,  my  slow  subsystem  will  be  (3.32)  with  n*  =  0.  In  order  to  derive  the 
fast  subsystem,  I  first  rewrite  the  error  system  (3.32),  (3.33)  in  terms  of 
a  fast  time  variable,  x  =  -^  t .  Then  I  set  p=0  which  corresponds  to  making 
my  nominal  plant  constant  at  its  initial  values.  The  resulting  fast  subsystem 
is  given  by 


’  A22"' 


(3.41) 


In  order  to  study  the  slow  subsystem,  I  form  a  system  of 

differential  equations  from  (3.32)  by  setting  n '  =  0  and  e  =  0 ,  and  then 
T  1 '  T  2  '  T  T 

replacing  [x'  ,v  ,v“  ]  with  X  and  8'  with  <f> .  This  results  in 


where 


A  = 


X  =  AX  + 

-  *  T  3n-2 

b (w  +  DX)  i>,  xe  ]Rjn  *■  , 

(3.42a) 

•k  -T 

?n 

b  =  -  r 

(w  +  DX)h  X,  4>*=  ®  , 

(3.42b) 

i 

i  *  T 

an  +  doblh1 

i  i 

!  bic  !  bid 

i  i 

*  T 

d  bh. 

0  1 

1  *T  1  *T 

j  A  +  be  j  db 

i  i 

(3.43) 

bh^ 

;  o  i  a 

i  i 

‘ 

-T  T  T 

b  =  [bj,b  ,0] 

» 

(3.44.) 

hT  =  [hi, 0,0] 

9 

(3.43) 

lx(3n-2) 

^(n-l)xn 

1 

1 

1 

1 

1 

Vi 

1 

1 

1 

1 

1 

^(n-l)xn-l 

D  = 

*T 

1 

j 

1 

1 

hT 

1 

1 

1 

__L_ 

°lxn-l 

1 

1 

1 

1 

_ _i._ 

°lx (n-1) 

^(n-l)xn 

1 

1 

1 

1 

1 

^(n-l)x(n- 

■u| 

Xn-1 

(3.46) 


Note  that  DX  replaces  w'.  From  Figure  3.3,  it  follows  that  the  transfer 
function  of  the  linear  time  invariant  part  of  (3.41)  is 


-T  -1- 

h  (si-  A)  b 


*  Z  (s)Z  (s) 

c  k  - - - 1 j - - - 

°  P  Z  (s) R  (s)  -  c*TT(s)R  (s)  -  k  (d  Z  (s) + d*TT(s))Z  (s) 
m  p  p  p  o  m  p 

(3.47) 


Theorem  3.8:  Let  r(t)  be  a  signal  which  satisfies  (2.17).  Let  u ^  =  miniu^y^} 
Then,  Vy  6  [0^2),  the  zero  solution  of  (3.42)  is  exponentially  stable  in  the 
large. 

Proof :  By  Theorem  2.2  the  transfer  function  (3.47)  is  strictly  positive  real. 
Then,  by  the  Kalman-Yacubovich  lemma,  a  positive  definite  symmetric  matrix  P 
exists  such  that 


ATP  +  PA  =  -  qqT  -  €  L  ,  L  =  LT>0,  £>0  (3.48) 

Pb  =  h  ,  (3.49) 


and  the  Lyapunov  function 


37 


has  Che  derivative  for  (3.42) 

W,,  .,.(X,$)  =  -  XT(qqT  +  £L)X<  0  .  (3 

(3.42) 

From  (3.50),  (3.51)  and  Theorem  3.2  it  follows  that  the  zero  solution  of 
(3.42)  is  uniform-stable. 

Because  <0,  every  solution  of  (3.42)  initially 

in  the  set  G  =  {X,b:  W(X,4>)  <  £}  remains  in  G.  Furthermore,  by  Theorem  3.5 


X->  0  as  t  -*• 00 


(3.52) 


Since  the  elements  of  w  are  bounded  continuous  functions  of 
time  with  bounded  derivatives  and  the  states  of  the  system  (3.42)  are  bounded, 
the  states  will  be  continuous  functions  of  time  with  bounded  derivatives. 
Hence,  the  elements  of 

C  =  w*  +  DX  (3.53: 

will  be  bounded  continuous  functions  of  time  with  bounded  derivatives.  Yuan 
and  Wonham  (1977)  have  shown  that  since  the  elements  of  :  are  bounded  con¬ 
tinuous  functions  of  time  with  bounded  derivatives,  (3.42)  and  (3.52)  imply 


C  b  0  as  t 


(3.54) 


From  (3.42b)  and  (3.52),  p  converges  to  a  constant  as  t-1-*. 
Yuan  and  Wonham  also  showed  that  if,  for  any 


aT;(t)  =0  t>0 


a  =  0 


(3.55 
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TO 

II 


ad 


then  <)>  ■*  0  as  Because  of  (3.52),  (3.55)  can  be  replaced  by  the 


condition 


a  w  (t)  =  0  t>0=>a  =  0 


A  sufficient  condition  for  (3.56)  to  hold  is  that 


X  * 

aw(jv)=0  -oo  <  v  <  ^  a  =  0 


&  5V 

where  w  (jv)  is  the  Fourier  transform  of  w  (t) .  Because 


m 


(3.56) 


(3.57) 


w  (jv)  = 


r(iv) _ 

kpRm(jv)Zp(jv)Wf(jv)  ’ 


(3.58) 


T(jv)Rp(jv) 
kpZm(jv)Zp(jv)Wf  (jv) 
kpT(jv)Zp(jv)Wf(jv) 


Theorem  2.1,  guarantees  that  condition  (3.57)  is  satisfied. 


Since  the  zero  solution  of  (3.42)  is  uniform-stable  and  every 

solution  of  (3.42)  with  bounded  initial  conditions  converges  to  zero  as  t-*-°°, 

the  zero  solution  is  unif orm-asymptotically  stable.  If  (3.42)  is  rewritten 
* 

with  c,  replacing  w  +  DX,  it  has  the  appearance  of  a  linear  system.  Hence, 
every  solution  of  (3.42)  with  bounded  initial  conditions  can  be  generated  by 
a  uniform-asymptoticallv  stable  linear  system.  Then,  by  Theorem  3.1,  each 
solution  of  (3.42)  which  begins  with  bounded  initial  conditions  decays 
exponentially  fast  to  the  zero  solution  and  the  proof  is  complete. 


Let  Y  =  [X  ,6  ]  .  Then,  from  Theorem  3.8  and  Definition  3.3, 


there  exists  and  >  0,  and  for  any  8 > 0,  there  exists  a  K^(S)  >  1  such  that 

if  Y(t;  Y  ,t  )  is  a  solution  of  (3.42),  then 
'  ’  o  o 

llY(t=  Yo-to)I:Ki(e)e'“1<t't°)!YoI  •  <3- 

By  Theorem  3.7,  there  exists  a  Lyapunov  function  V^(t,Y)  which  has  the 
following  properties: 

|]Y||  <V1(t,Y):K1(3)||Y||  ,  |Y|<3  ,  (3. 

|v]_(t,Y1)  -V1(t,Y2)|  <  L(t,S)| Y1-Y2|  ,  || Y1||  <  B,|Y2||  <  0  ,  (3. 

D0  42)Vl(t,Y)  q  aiVl^’Y)’  where  0  <  q  <  1  .  (3. 

Thus,  the  existence  of  a  Lyapunov  function  for  the  slow  subsystem  is 
established . 

Define  a  such  that  the  real  part  of  each  eigenvalue  of  A22 
is  less  than  or  equal  to  -a,  that  is 


a  =  inf{-Re[A(A22) ] } 


(3. 


It  follows  from  (3.63)  and  the 
there  exists  a  K0  such  that  if 


linear  time 


n’(x;  n^,To) 


invariant  nature  of  (3.41) 
is  a  solution  of  (3.41), 


that 

then 
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Replacing  x  with  in  (3.64)  and  applying  Theorem  3.6  shows  that  there  exists 
a  Lyapunov  function  V2(t,n')  which  has  the  following  properties: 

| n ' |  <  v2(t,n’)  <  k2| n’l  ,  (3.65) 

|v2(t,n”)  -  V2(t,n')|<  Kjn'  -  n"  J  ,  (3.66) 

°"(3.54)V2(t’n')  -  ~  v  v2(t’n,)  •  (3.67) 


Thus,  a  Lyapunov  function  for  the  fast  subsystem  has  been  found. 

Before  composing  these  two  Lyapunov  functions  to  form  a 
Lyapunov  function  for  the  error  system  (3.32),  (3.33),  I  need  to  simplify 
the  expression  for  u'.  Let 


X  = 


[x 


T  1,T  2' 

,v  ,v 


T 

J 


T 


(3.68) 


Y  =  [X 


T  <  T  T 

A  I  Q  .  1  1 


(3.69) 


With  e  =0,  (3.37)  can  be  rewritten 


•  T  *_T  A  _T  T  T  _T 

u'  =  -  w*  Tw  h  X  -  2w  TDXh  X  -  X  D  TDXh  X 


aT  j  ATT  T  31 

+  w  9'  +  f.Y  +  d  hXhV  +  Y  9' 

1  o  1  1  2 

+  d'h^b  h  n'  +  d  ' hub,  9  ' "^w' 
o  1  1  2  oil 


(3.70) 


Because  r  and  r  are  bounded,  the  elements  of  w  and  w  are  bounded  and 
there  exist  constants  G^,  ,G2  ,G^  ,G^  ,G^  €  such  that 


1  -  -L  ,  1GilYl  +SlYl  +G3lYi  +  GJn'I  +  G sHM1  •  (3-71) 

K2IA22b2l 


The  main  result  of  this  section  is  contained  in  Theorem  3.9. 

Theorem  3.9:  Let  r(t)  be  the  same  r(t)  as  in  Theorem  3.8.  Then,  there  exists 

4-  k 

a  such  that,  Vu  €  [0,^)  the  zero  solution  of  (3 . 32) ,  (3 . 33)  with  e  =  0  is 

exponentially  stable. 

Proof :  Take  u  ^  S  u  2  *  ^en’  t^ie  Lyapunov  function  V^Ct.Y)  exists  and  has 

properties  (3.60)— (3-62) .  Using  these  properties  and  (3 . 65) -(3 . 67) ,  I  can 


write 


D(3.32)e*=0Vl(t’Y)  -  “  q  aivi(t»Y)  +  K3(6)V2(t,n')  ,  ||Y||  <  6  ,  (3.72) 

Dn.33)e*-0V2(t’n’}  (7  "  G4  "  G5V1  (t  ,Y)  )V2  (t  ,n  ' )  +  f  (Vj.  (t  ,Y) )  t  ,Y) , 

(3.73) 


where 


K3(S)  =  sup{L1(t,6)|b1J|h2| 


f(V1(t,Y))=  G±  +  G2V1(t,Y)  +  G3V^(t,Y) 


Define  M(Y,n’)  as  the  solution  of 


(3.74) 


(3.75) 


M(Y  ,  n  ' )  =  K. 


uK3(M(Y,n’)) 

Yl  +  a  -  u (G4  +  q  a1  +  G,M(Y , n ' ) )  K2l n  I 


(3.7b) 


which  satisfies  M(Y,n  '  )  =  (||  Y|j  )  ||  Y||  +  ^  K3(M(Y,n '  )  )K„ 


Let  M  =  M(Y  ,n ' )  and  define 
o  00 

uK3(Mq) 

a_  =  q  01 - - - — - -  f  (M  ) 

3  1  a  -  u  G,  -p  GCM  -u  q  a,  o 

4  5  o  1 


(3.77) 


^ _ jl_ 


Lm.  i-.  _  _ 


Next,  consider  the  Lyapunov  function 


PK  (M  ) 

V,(t,Y,n')  =  V.(t,Y)  + - - - — - 

3  1  a  -  p  G,  -  p  G^Mq  -  p  q  a1  2 


V0(t,n')  • 


(3.78) 


The  right-side  derivative  of  V^(t,Y,n')  along  the  trajectories  of  (3 . 32) -(3 . 33) 


with  e  =  0  satisfies 


*(3.32)-(3.33)e*  H  0V3(t,Y,n,)  -  “[q  “l  o- 
o-P  G4-  pG5V1(t,Y) 


pK3(mo) 


p  G4  -  pG5Mo  -  pq 


•  f(V1(t,Y))]V1(t,Y) 


3  v“o'  va  -  p  G.-pGcM  -  pq  a.  “1>v2  (t ’ n ’  > 
4  3  0  1 


-  K,(M„)(- 


=  -  a3V3(t,Y,n’) 


PK  (M  ) 

„  -„r  M - — — (f(M  )  -  f(V,(t,Y))V  (t,Y) 

a-y  G.  -  yGcM  -  yq  a-  o  I  1 

4  _>  o  1 


PK  (M  )G 

o-p  G/  -  pG.M  'I  pq  «V(Mo  -  V1(t,Y))V2(t,n') 


4  5  o 


-  c 


PK3(mo) 


a-p  G.  -  pGcM  -pq  a.'  "“o7'? 
4  5  o  1 


rf(M„)v,(t,n’)  . 


From  (3.60),  (3.65),  (3.76)  and  (3.78)  it  follows  that 


±  Mo 


From  (3.78)  it  is  clear  that 


v1(t,Y)  <  v3(t,Y,n') 


From  (3.79),  (3.80)  and  (3.81),  it  follows  that  if  ^3 >  0,  then 


a,t 


V,(t,Y,n')  1  V_ (0,Y  ,n  '  )e_  3 


(3.79) 


(3.80) 


(3.81) 


3 


(3.82) 


Hence,  defining 


oq  ot1 

U4  =  (q  aL)  (G4  +  q  o^)  +  K3(0)G1  (3’83 

and  taking  =  minlv^l1^}  completes  the  proof. 

An  estimate  of  the  region  of  attraction  is  available  from  the 

k  -f- 

proof  of  Theorem  3.9.  For  every  [O,^)  there  exists  Mq(u)£]R  such  that 

k  k 

if  Mq  S  [0,M  (u)),  then  a^>0,  and  if  Mq  =  Mq(m),  then  a^=0.  Then  the  region 
of  attraction  contains  the  set  D(p)  defined  as 


D(u)  =  { (Y,n)  :  M(Y , n ' )  <  M*(p)} 


(3.84) 


From  (3.77)  ^it  is  clear  that  if  K^(6)  =  0(1)  for  all  B€1R+,  then  as  p-+0 

*  ~  + 
Mq(ij)->0(u  )  •  If  one  further  assumes  that  K^(8)  =  0(1)  for  all  B  e  1R  , 

then  as  y  -*■  0  the  region  of  attraction  includes  points  which  satisfy 


IKl  ; 

II  noI  : 0<M 


Thus,  the  zero  solution  of  (3 . 32) -(3. 33)  with  e  =0  and  u  ^small  enough  is 

~o* 

exponentially  stable  and  the  region  of  attraction  is  0(u  )  in  the  slow 

.2 

variables  and  0(y  ^ )  in  the  fast  variables. 

Remark  3.2:  In  this  section  I  have  required  that  r(t)  have  exactly  n  distinct 

sinusoids  in  the  low  frequency  range.  However,  the  "exactly  n"  part  of  this 

k  k 

requirement  was  used  only  to  choose  9  (y) .  If  we  choose  9  (u)  by  some  cri¬ 
terion  other  than  making  the  steady  state  value  of  the  tuned  error  equal  to 
zero,  then  we  can  replace  "exactly  n"  with  "at  least  n"  and  the  analysis  of 
this  section  will  still  hold. 


3.4.  Effects  of  e  #  0 

* 

In  this  section  I  show  that  the  effects  of  e  r  0  are  only 

0(u) .  This  is  accomplished  in  two  steps.  First,  I  show  that  I  can  choose 

initial  conditions  on  the  tuned  system  so  that  |e  ]  <  0(p)  for  all  t >  0. 

* 

Next  I  show  that  if  Mo  represents  the  boundary  of  the  region  of  attraction 

■k  & 

when  e  =0,  then  the  boundary  of  the  region  attraction  when  j  e  |  =  0(p)  will 
be  represented  by  (1-0  (uz  ))Mq.  At  the  same  time  I  show  that  solutions  of 

i  *  i 

(3. 32)— (3. 33)  with  | e  |  —  0(u)  beginning  in  this  slightly  reduced  region  of 
attraction  converge  to  an  0(y)  residual  set  containing  zero. 

In  order  to  show  that  I  can  choose  initial  conditions  on  the 
tuned  system  so  that  |e  |  <  0(p)  for  t  >  0,  I  begin  by  choosing  a  different 
representation  of  the  model.  Define 
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*  T 

+  0  (0)  D  . 

(3.85) 
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i 
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Then  the  model  equations  (3.20)— (3.21)  can  be  replaced  by 


X  =  AX  +  be  (0) r  ,  X  (0)  =  X 
m  m  o  m  m 


-T 

v  =  h  X 
'  m  m 


Next,  I  find  a  Lyapunov  function  for  this  representation. 


(3.86) 

(3.87) 


Since 


T  .  Z  (s)Z  (s)Z(s) 

hTrsT  _  A)~Lhr  =  k  — _ E _ 2 _ 

h  (si  A)  bcQ  m  z  (s)z  (s)R  (s) 

m  p  m 


and  because  Z  (s) ,  Z  (s)  and  R  (s)  are  Hurwitz,  I  know  that  A  is  a  stable 
m  p  m 

matrix.  Define  so  that  the  real  part  of  each  eigenvalue  of  A  is  less 
than  -a  , 

a1  =  inf {-  Ref  X ( A) ]  }  .  C 

Then,  if  X  (t;  X  ,t  )  is  a  solution  of  (3.86)  with  r  =0,  there  exists 
m  mQ  o  ’ 

a  K  such  that 
4 


X  (t;X  ,t  ) 
m  m0  o 


K  e"al(t  to)||  x  |i  ,  t  >  t 
4  <  “  o 


By  Theorem  3.6  and  (3.89),  there  exists  a  Lyapunov  function  V  (t,X  )  such 

4  m 


i  VC'V  :  XJ 


D(3.86)rE0V4(t-V  i  -  °lV4(t’V 


The  first  result  of  this  section  is  contained  in  Theorem  3.10 

*  *T  i*T  2*T  T 

Theorem  3.10:  Let  X  =  [x  ,v  ,v  ]  .  Require  r(t)  and  r(t)  to  be 
*  *  * 

bounded.  If  X  (0)  =  X  (0),  n  (0)  =  0,  and  uS[0,u,),  then  there  exists 

m  3 

ill,  €  jR+  such  that  for  t  >  0 


Then  the  Lyapunov  tunction 


*  u  III, 

v5(t,E,n)  ■  V4U, E)  + 


v  (t  ,n  ) 

uo,  2 


(3.103) 


satisfies 


(3 


.93)(3.94)V5(t’E’n  }  -  ‘  ^V5(t,E,-*)  +  u(K4m5ni2  + 


K4E7K9an1  +  m2m3) 


o  -  pG. 


1 


(3.104) 


where 


k.iilk-iil 

T  .  ij  ( K  m .  +  -2t~g~g . 3  - 

1  4  ^  a  -  uG.  -  pa. 


(3.105) 


If  h€[0,u3),  then  a.  >  0.  Hence  E*(0)  =  0,  n*(0)  =  0,  (3.90),  and  (3.104) 


imply 


|e)  ■'  V,  ( t ,  E)  -  V  (t,E,n  ) 


3 


k,il,  K0(nu  +  to^rl) 
(m  +  i- - ^  ) 

Xr  ■>  0  -  uG.  -  UO, 


(3.106) 


The  proof  is  completed  by  defining 


\  ■  I  hi 


K,m0  K.jdn.  +  m7m9) 

(ms  +  -r—k — ) 


3  o  -  uG.  -  uo, 

4  1 


(3.107) 


Now  that  I  have  a  bound  on  !  e  |  ,  I  examine  the  effect  of  this 
bounded  input  on  solutions  of  (3. 32)-(3. 33) .  The  approach  is  similar  to  that 
when  e  =  0.  First,  I  bound  iu'|  .  Then,  I  write  inequalities  for  d"!" 


VL(t,V)  and  ^ 


(3.32) 

(t ,n ' ) .  Next,  I  define  some  terms  needed  to  state  the 


final  result  of  this  chapter.  Finally,  I  summarize  the  results  of  this 
chapter  in  Theorem  3.11. 

As  stated,  I  first  find  a  bound  on  [u'j  assuming  that  e 


Because  it  w  i  is  bounded  for  all  time,  there  exist  constants  G.  and  0  such 


48 


_  ju 

K2ll A22b 2 1 W  rw”  1  g6  »  t  >  0  »  (3.108) 

G1  +  p2K2||  A~2b2|[  ||  rw*J B6  <  Gx  ,  t  >  0  .  (3.109) 

Next,  define 


G2  G2  +  !jK2lA22b2ll  N®€ 


(3.110) 


Then,  I  can  bound  Ju'j  for  t  >  0  by 


u  :  < 


K2| A22b2 


rrr  [5illYl  +  SH2  +  S!v!3  +  G4IMI  +  SUIKI 


yc6ni6] 


(3.111) 


Now,  I  need  to  write  inequalities  for  D(3.32)Vl(t,Y)  and 
0^  33)  -and  define  terms  I  use  to  state  the  final  theorem.  Using 

(3.60)-(3.62) ,  (3 .65)-(3 .67) ,  and  (3.111)  I  write 


D0.32)Vl(t,Y)  <  -  a2(8)V1(t,Y)  +  K3 (g)V2 ( t , n  ’ )  +  uK  (8)Bg  , 

Dp.  33)V2(t»r>’ )  <  -  G4-G5V1(t,Y))V2(t,n*)  +  f  (V1(t,Y))V]_(t,Y) 

+  mg6ih6 

where 


(3.112) 


(3.113) 


o7(S)  =  qal  -  uB6|rjK3(B) 


(3.114) 

(3.115) 

(3.116) 

Recall  that  for  the  e  =  0  case  I  defined  M(Y,n')  which  was  used  to  bound 


K.(g)  =  sup  J'L1(t,8)||  Fw  ( t )  ||  } 
3  t>0 


f(V1(t,Y))=  G1  +  G2Vx(t,Y)  +  G;jV1(t,Y) 


•I Y ||  and  estimate  the  region  of  attraction.  In  this  case  M(Y,rv)  will  not 
always  be  a  bound  on  ||y||,  and  therefore,  I  define  as  the  solution  of 


Mi  = 


K  (pM  )G 

fll,  [K  (pM  )  +  p  - ,  =-  „ - T] 

6l  5V  1  o  -  p[G^  +  pG^M^  +  ] J 


K3(uM1)f(pM1) 

CT2^:1l)  a  -  p [G^  +  uG^M^  +  qa^] 


(3.117) 


.  a6V0) 

which  satisfies  M,  =  - r— ; — 

1  ao(0) 


Next,  I  define 


M2  =  max{M(YQ,n^) ,pM1) 


(3.118) 


which  I  will  use  to  bound  |Y| .  The  rate  of  exponential  decay  will  be  greater 
than 


K3(M2)f(M2) 

*6  =  °2(M2)  ‘  U  a  -  p[Ga  +  G5M2  +  qn,  ] 


(3.119) 


Define  u,  as  the  solution  of 
4 


0  =  °2(0)  '  M4  7 


k3(0)g1 


“4[G4  +  9“!^ 


(3.120) 


which  satisfies  =  p^,  where  p^  is  defined  in  the  proof  of  Theorem  3.9. 

*  —  4- 

To  estimate  the  region  of  attraction,  define  a  function  M0  :  [ 0 , u  ,  ]  *>  IR  such 

l  4 

*  : k 

that  if  M2  <=  [0,M2(p)),  then  > 0  and  if  M2  =  M?(;j),  then  <  0.  Next, 


define 


(3.121) 


K3(M2)G6 


M3  <*6  [K5(M2)  +  u  a  -  u[G4  +  G5M2  +  qc^]1 


GA 

M  =  - - — - - 

4  a  -  p[G,  +  G5M2] 


f(pM3)M3 


a  -  u[G4  +  G5M2  +  3qct3  ] 


(3.122) 


Finally,  I  define  three  sets  which  are  useful  for  describing  the  behavior  of 
solutions  of  (3.32)-(3.33) , 

D1(u)  =  {Y  ,n '  :  M(Y.n')  <  mM3,  Kjn'H  <  U2M4)  ,  (3.123) 

D2(p)  =  { Y , n ’  :  |  Yj|  <  uM3 , 1  n  ’ [(  <  ,  (3.124) 

D3(u)  =  (Y,n’ :  M(Y ,n 1 )  <  M*(p)}  .  (3.125) 

The  set  D3(u)  will  be  an  estimate  of  the  region  of  attraction,  the  set 
D2(u)  will  be  a  residual  set  to  which  all  solutions,  beginning  in  D3(y) , 
converge  exponentially,  and  the  set  D^(p)  is  a  region  such  that  solutions 
beginning  in  ^^(p)  stay  in  D2(u)  for  all  t  >  0.  Theorem  3.11  shows  that,  for 
y  small  enough,  these  statements  are  true.  Theorem  3.11  also  provides  a  more 
detailed  bound  on  JY|  and  j|n'||  . 

Theorem  3,11:  Let  r(t)  be  the  same  as  in  Theorem  3.8.  Let  (Y  ,n ' ) (t ;Yq, n^) 
be  a  solution  of  (3 . 32) - (3 . 33)  beginning  at  (YQ,n^)  at  time  t  =  0.  Then, 
there  exists  a  P  €  ]R+  such  that  for  all  uG  [0,“),  if  (Yq,V)£  ®3(u)  »  then 


I Y II  <  M(Y  ,n')e  a6c  +  uM_  (1  -  e_a6C) 
II  ii  -  o  o  3 


(3.126) 


ni  :  K2  noi  exP[-<u  -  G4  -  SM2)t] 


o  -  y  [G.  +  G,-M„  +  qa1  ] 


[1  -  exp[-(-  -  G,  -  G,M0)t] 


(3.127) 


f(M(Yo,n;)e"a6c  +  um3(1  -  e~a6t))  [M(Yo,n;)e  ^d-yM^l-e  a6C)] 


a  -  y  [G^  +  G  M2  +  Iqo^] 


for  t  >  0 . 


Proof :  The  proof  takes  several  steps.  First,  I  define  y  so  that  for  all 
y€  [0,  y)  D^(y)CD^(y).  Next,  I  assume  a  bound  on  JjYjj  exists  and  show  that 
if  (Yq,V)  €  D3(u)  ,  then  M2  is  a  bound  on  jj  Y||  .  Confirmation  of  M2  as  a 
bound  on  ||  Y||  leads  directly  to  the  bound  (3.126).  Finally,  I  use  the  bound  on 
II Yll  to  get  the  bound  (3.127)  on  ||n'||. 


Because  K^(8)  is  a  nondecreasing  function  of  B,  it  follows  that 

*  —  k 

M2(y)  is  a  decreasing  function  of  y  on  y€~  [0,y^].  Furthermore,  M2(0)  =  +  » 
and  M2(y^)  =  0.  Clearly,  there  exists  a  (0,u^)  such  that  if  u£  [O.y^), 

k  k  — 

then  M^Cu)  >  and  ^2^5^  -  U5M1#  Ta^e  V  ~  mini  } .  Now, 

take  yG  [0,  y),  and  assume  that  a  constant  exists  which  bounds  |Y|  for 

t  >  0.  Define 


k3(m5) 


VsU.Y.n')  .  V^t.Y.n')  +  »  —  P|  +  '  -  V2(t,n')  .  (3.128) 


Using  (3.111)  and  (3.112)  I  can  write 


uK3(M5)f(M5) 


D(3.32)(3.33)V6(t,Y,n,)  -  “(a2(V  “  a-u[G4  +  G^  +  qa1])Vl(t,Y) 


mK3(m5) 


a  -  u[G4  +  G5M5  +  qct1  J  qalV2 (t ,n ' > 


uK  (M  )G 

+  uflUK.  (M_)  +  - — -j— 2  P—  -  i 

6  4'  5  a-  y[G4  +  G5M2  +  qc^] 


(3.129) 


<  -  (a2(M5)  -  u 


K3(M5)f(M5) 


a-  p[G4  +  G5M5  +  qa1])V6(t,Y,n’) 


UK  (M  )G 

+  “WM5>  +  ~^G4V|-5M-3---rqai]^ 


Assume  further  that  <  M2(u).  Then 

K  (M  )f(M  ) 

V6(t,Y,n’)  <  V6(t.Yo,n;)  exp[-(a2(M5)  -  U  qa,]^ 

K  (M  )G 

VW  +p  a  -  u[G4  +  G5 M5  +  q 
+  U  IT  K3(M5)f(M5)  ‘ 


(3.130) 


CT2(M5}  P  a-  p[G4  +  G5M5  +  qc^] 


In  order  for  to  bound  |Y|  I  can  choose  to  satisfy 
► 

UK  (M  )K  In'! 

M.  =  max  K.(|Y  |)||Y  I  +  - j— — — ~ -  , 

D  )  1  II  oil'll  oi  a-u[G4  +  G^M^  +  qa^] 

K-(M,)G, 

K3(M5)f(M5) 


(3.131) 


°2(M5)  '  U  a  -  u[G.  +  G_M_  +  qa,  ] 


If  Mg  is  equal  to  the  first  term  on  the  right-hand  side  of  (3.131),  then 

Mg  =  MCY^n^).  If  is  equal  to  the  second  term  on  the  right-hand  side 

of  (3.131),  then  Mg  =  yM^.  Hence,  Mg  =  M2  and  M2  bounds  j|Y|  for  t  >  0. 

£ 

For  (Y  ,n^) €  Dg(u) ,  then  M2  <  M2(u)  and  I  can  write 

| Y||  <  V1(t,Y)  <  V6(t,Y,n')  <  M(Yo,n(;)e“a6t  +  yM3(l-e“a6t)  (3.132) 

which  provides  the  bound  (3.126)  on  ||Y||. 

Next  consider  the  function 


uf(V  (t,Y))V  (t,Y) 

Vy(t ,Y,n  )  =  v2(t,n')  -  c  _  p[G^  +  g5m2  +  3qot  1  ] 


(3.133) 


Using  (3.112)  and  (3.113),  I  write 


D(3.32)(3.33)V7(t’Y^’>  -  “  °4  ‘  SW  C  ’ n  ’  >  +  ‘Vfi 


ya9(M?)(G.  +  2G  V  (t,Y) +  3G  V^(t,Y)) 

+  +  a,2u[^  +\--M— ^ - 1  v, ( t »y) 


-  -  (y  ’  G4  -  G5M2)V2(t*n,)  +  uG6m6 


(3.134) 


u3c  (M  ) 

+  +  q-M[G4-bG5M2^3qa111Vl(t>Y) 


-  (-  -  G,  -  G.M„) V, ( t ,Y , n ' )  +  uGcHL 


From  (3.134),  it  follows  that 


V7(t,Y,n’)  <  V7(0,Yo,n;)  exp[-  (2.  -  G4  -  G^Jt] 


a  -  u[G4  +  G5M2  +  3qa;[] 


[1-  exp[-  (-  -  G4  -  G5M2)t]] 


Using  (3.65),  (3.132),  (3.133)  and  (3.135)  leads  to  the  bound  (3.127)  and 
completes  the  proof. 

In  the  introduction  to  this  section,  I  claimed  that  the  effect 
of  the  0(p)  e*  on  the  region  of  attraction  was  very  small.  To  see  this, 
assume  that  K  (8)  =0(1)  for  all  6€R+,  || T ||  =0(1),  and  sup|u*(t)  ||  -  0(1)  . 

J  t 

Under  these  assumptions,  M*(u)  is  chosen  so  that  if  Mq£;  [0,M*(u))»  then 


VV 


0  <  qa  -  -  u  - - rr  j.r  v j.  i  £(S 

1  7  o  O  -  u[G4  +  G^M^  +  qa^J  o 


where  M  >  0  is  an  0(1)  valued  function  of  (M  ) .  Recall  that  M*(u)  was 
/  o  o 


chosen  so  that  if  M  €(0,M*(u))>  then 

o  o 


K,(M  ) 

°  qai  '  U  a-u[G  +G  M  +  qa  ]  f  (Mo*  ‘ 
4  5  o  1 


Comparing  (3.136)  and  (3.137)  we  see  that  as  p  +  0 


M*(u)  =  (M*(u)  -  0(u) ) 

z  o 


and  thus  we  have  M*(u)  = M*(u) (1  -  0(u  ))• 
z  o 

Remark  3.3:  Notice  that  Theorem  3.10  does  not  restrict  the  number  of 

sinusoids  in  the  reference  input.  That  is,  the  bound  j  e*  I  <  ufll  does  not 

0 

depend  on  the  number  of  sinusoids.  Using  a  relaxed  version  of  Theorem  3.8, 
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as  suggested  in  Remark  3.2,  one  could  relax  the  hypothesis  of  Theorem  3.11 
to  r(t)  sufficiently  rich  and  low  frequency  without  changing  the  conclu¬ 
sions  of  Theorem  3.11.  If  we  relax  the  "exactly  n”  requirement  then  we 
must  choose  6*(p)  by  some  criterion  other  than  making  the  steady  state 
e*  =  0.  An  interesting  criterion  is  9*(u)  =6*(0).  For  this  choice  of  9*, 
the  previously  discussed  relaxed  version  of  Theorem  3.11  would  be  an  excellent 
proof  that  low  frequency  sufficient  richness,  i.e.,  dominant  richness, 
provides  robustness  with  respect  to  high  frequency  unmodelled  dynamics. 

I  use  robustness  here  in  the  sense  that  in  the  presence  of  high  frequency 
unmodelled  dynamics,  the  adaptive  control  system  remains  stable  and  tracking 
and  parameter  errors  become  and  remain  small,  0(p). 


3.5.  Discussion 

In  this  chapter  I  have  shown  that  trajectories  of  the  model 

reference  adaptive  control  system  excited  by  n  low  frequency  sinusoids 

which  begin  in  a  region  of  attraction  will  converge  exponentially  to  an  0(u) 

residual  set  if  u  is  sufficiently  small.  Furthermore,  if  r(t)=r*(t),  the 

steady  state  value  of  e*  is  zero  and  the  equilibrium  eventually  becomes 

exponentially  stable.  It  was  also  shown  that  if  the  rate  of  convergence  of 

the  unperturbed  'IRAS,  i.e.,  when  u = 0,  is  0(1),  then  the  region  of  attrac- 
_  1  _  3 

y  y 

tion  may  be  0(u  L  )  in  the  slow  variable  and  0(u  )  in  the  fast  variables. 

Through  Remarks  3.2-3. 3  I  have  also  detailed  a  method  by  which  one 
can  prove  that  robustness  with  respect  to  high  frequency  unmodelled  dynamics 
can  be  gained  via  low  frequency  sufficient  richness  of  the  reference  input. 
This  method  was  based  upon  the  idea  of  fixing  the  reference  input  and  then 
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showing  that  for  sufficiently  small  y,  exponential  stability  is  retained. 

How  does  one  use  the  ideas  in  this  chapter  if  y  is  fixed  and  the  problem 
is  to  choose  r(t)  and  T  to  provide  robustness  with  respect  to  the 
perturbation?  The  first  step  is  that  a ^  must  be  positive.  From  (3.75) 
and  (3.77)  we  see  that  at  the  very  least  one  needs  <  0(y  *) .  This  leads 
for  one  condition  to  w  rw*<0(y  ).  Another  requirement  might  be  that 

_  i  _  I  _  I 

m6<0(y  2)  .  This  leads  to  |r|  <  0(y  ^)  and  |r|  <0(y  ^)  .  We  begin  to  see 

that  keeping  ||r|,  |r|,  and  |r|  all  small  with  respect  to  y  2  wm  probably 

accomplish  the  desired  goals. 

Before  going  on  to  examples  in  Chapter  4,  I  would  like  to  discuss 
the  relationship  of  this  work  to  that  of  other  authors. 

Remark  3.4:  While  Remark  2.3  suggested  that  dominant  richness  was  important 
for  the  existence  of  a  unique  equilibrium,  it  is  easily  seen  that  the  sta¬ 
bility  of  this  unique  equilibrium  requires  a  stronger  dominant  richness 
condition.  For  a  given  y > 0,  each  w^  in  the  input  must  be  less  than  0(y  *) 
to  assure  existence.  As  just  discussed,  stability  can  be  assured  only  when 
the  existence  condition  is  met  and  when  each  product,  r^w^,  is  less  than 
0(y  ).  This  more  restrictive  condition  corresponds  exactly  with  the  notion 

of  dominant  richness  as  discussed  by  Ioannou  and  Kokotovic  (1983)  and  Krause 
(1983). 

Remark  3.5:  Krause  (1983)  has  also  found  that  when  the  reference  input  is 
dominantly  rich  and  y  is  sufficiently  small,  the  convergence  of  the  trajec¬ 
tories  of  the  model  reference  adaptive  control  system  toward  e  =  0,  6  =  9*(0) 
for  0(1)  parameter  errors  can  be  assured  by  choosing  ||r||  sufficiently  small. 
The  results  of  this  thesis  agree  in  the  following  way.  If  the  assumptions 


of  dominant  richness  and  y  sufficiently  small  guarantee  that  a,  > 0  when 

6 

r  =  0,  then  there  must  exist  T  with  |jr|J  small  enough  so  that  a  >0  with 

6 

II T ||  >  0 .  Since  >  0  is  sufficient  for  the  convergence  of  trajectories  in 
D^(u)  to  D  (u) ,  the  trajectories  will  then  behave  as  previously  described. 
Remark  3.6:  While  it  is  easy  to  see  that  the  exponential  stability  of  this 
system  with  e*  =  0  can  be  used  with  standard  stability  results  to  show  robust¬ 
ness  with  respect  to  additive  disturbances,  it  is  not  clear  what  happens  if 
the  reference  input  is  the  sum  of  a  dominantly  rich  input  and  a  high  frequency 
input.  Astrom  (1983)  showed  via  averaging  techniques  that  for  small  ||r|,  the 
low  frequency  excitation  provided  robustness  with  respect  to  the  high 
frequency  excitation.  This  can  be  partially  seen  by  realizing  that  e*  is  the 
difference  of  the  outputs  of  two  low  pass  filters,  and  hence,  |e*|  will  be 
very  small  for  components  of  the  reference  input  which  are  high  frequency. 

Thus  the  forcing  effect  of  high  frequency  inputs  for  slowly  varying  9(t)  will 
be  small.  However,  the  techniques  of  this  thesis  may  not  be  able  to  prove 
this  because  of  the  w*  coupling  term  between  fast  and  slow  subsystems. 

Remark  3.7:  As  a  final  comment  on  this  chapter,  I  must  say  that  the  results 
of  this  chapter  could  be  interpreted  as  a  special  case  of  the  work  done  by 
Kosut,  Johnson,  and  Anderson  (1983).  However,  I  believe  that  much  can  be 
gained  by  dealing  with  the  special  case  I  have  chosen  and  by  using  the 
Lyapunov  approach,  which  is  not  the  approach  of  Kosut,  et  al.  First,  the 
special  case  of  high  frequency  unmodelled  dynamics  which  1  have  chosen  is 
in  itself  a  fairly  general  and  very  important  topic  of  current  research  in 
adaptive  control.  Second,  the  representation  I  was  able  to  choose  for  the 
unmodeled  dynamics  allowed  for  a  natural  decomposition  of  the  problem 
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which  in  turn  provided  the  starting  point  for  the  composite  Lyapunov  function 
proof  of  the  results.  Finally,  the  results  help  to  point  out  that  high 
frequency  inputs  and  high  adaptation  gain  are  the  main  problems  to  watch 
out  for  when  applying  model  reference  adaptive  control  to  a  plant  with 
unmodelled  high  frequency  dynamics. 
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CHAPTER  4 
EXAMPLES 

In  Chapter  3,  I  derived  conditions  for  stability  of  a  model 
reference  adaptive  system  with  unmodeled  high  frequency  dynamics.  In  this 
chapter,  I  present  a  sample  system  of  this  type  and  apply  the  results  of 
Chapter  3  to  derive  conditions  for  stability  of  this  sample  system  with  two 
different  sets  of  parameters.  While  these  two  special  cases  will  not 
illustrate  every  aspect  of  the  theory  presented  in  this  thesis,  they  provide 
sufficient  insight  to  appreciate  the  method  suggested  by  the  theory  and  to 
discover  the  main  drawbacks  of  the  suggested  approach. 

This  chapter  is  divided  into  several  sections.  In  the  first 
section  I  present  the  sample  system  and  derive  the  associated  error  system 
following  the  ideas  of  Section  3.2.  The  second  section  is  devoted  to  a 
study  of  the  sample  system  with  the  parameters  chosen  so  that  it  becomes 
a  linear  time  invariant  system.  In  the  third  section,  the  parameters  are 
chosen  so  that  the  system  becomes  the  simplest  nonlinear  system  possible 
and  the  study  is  repeated  for  the  new  case.  In  the  final  section  I  discuss 
the  problems  illustrated  by  the  two  special  cases  and  some  aspects  of  the 
theory  which  are  not  illustrated  by  the  two  special  cases. 


4.1.  The  Sample  System 

In  this  section,  I  present  a  somewhat  general  model  reference 
adaptive  system  with  high  frequency  unmodeled  dynamics.  Then,  following 
Section  3.2,1  derive  the  error  system  representation  of  the  sample  system. 


This  error  system  is  used  with  two  different  sets  of  parameters  to  create 
the  examples  used  in  Sections  4.2  and  4.3. 


A  block  diagram  of  the  sample  system  is  shown  in  Fig.  4.1. 
The  differential  equations  describing  it  are 


model : 


plant : 


control: 


parameter 

update: 


x  =  -x  +  r , 
m  m 


x  =  ax  +  b 


45 
2  0 


z. 


-1 

.5 

0 

uz  = 

z  *f 

/5 

-.5 

-1 

/  2 

u  = 

fiT 

6  w. 

w 

u 

<111 

x]  , 

o  =  _ 


w(x-x  ) 
m 


c 

(4 

(4 


The  first  step  in  the  derivation  of  the  error  system  is  to  make  the 
transformation 


/I 

=  z  -  /  5 


n  =  z 


The  plant  equations  (4.2)  become 

x  =  ax  +  bu  +  b 


17!  «! 

n. 

(4.6) 

Next,  define  the  tuned  system 


.  * 

X 


(a  +  b6*)x* 


* 


+  b@^r  +  b 


where 


*  .*T  * 

u  =  0  w  . 


*  r 
w  =  [r 


.*iT 


0 


0 


The  tuned  error  is  then  defined  as 

*  * 


e  =  x  -  x  . 

m 


(■ 


Figure  4.2  shows  a  block  diagram  of  the  tuned  system  and  the  tuned  error 
Finally,  making  the  definitions 

i  * 

X  =  x-x  , 

e'  =  e-0* ,  ( 

i  * 

n  =  n-n  > 


I  write  the  error  system  equations 


=  ax '  +  bu '  +  b 


/!  " 


t 

n  , 


un  =  - 


h 


w(x 1  4-  e*)  , 


Figure  4.2.  Block  diagram  of  the  tuned  system  and  tuned  error  for 
the  sample  system. 
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This  corresponds  to  a  plant  whose  "slow"  pole  matches  the  model  pole  but 

which  has  an  unknown  dc  gain  of  y.  The  choice  e*  =  2  provides  the 

correct  dc  gain.  The  choice  of  the  adaptive  gain  =  1  is  arbitrary  and, 

hopefully,  reasonable  for  this  problem. 

One  of  the  first  things  to  notice  about  this  system  is  that  u* = 2 

is  a  constant.  Therefore  u* = 0  and  hence  it  follows  from  (4.1),  (4.8),  and 

(4.10)  that  the  choice  x*(0)  =x  (0),  n*(0)  =0  results  in  e*  =  0.  Using  this 

m 

and  (4.17)  I  write  u'  for  this  case, 


•  t 

u 


[x 


(4.18) 


i 

■ 


L 

> 


Rewriting  the  error  system  using  the  appropriate  modifications  for  this 

'  i  T 

special  case  and  Y =  [x  0^]  results  in 


(4.19a) 


■5-rii 

-1 

Ui 

_ _ _ 1 

pn  =  u  ; 

'  5  1  2 

[1  -.5]Y  + 

.'•5 

-i ! 

(4.19b) 


Now  I  want  to  apply  the  method  developed  in  Theorem  3.9  to  determine 
a  u  such  that  (4.19)  is  exponentially  stable  for  all  p£[0,u).  The  first 
step  is  to  consider  the  subsystems 


Y  = 


-1 


•5l 

o ! 


-i 


(4.20) 


Solutions  of  (4.20)  satisfy 


||y||  <  2 . 44e  ’  5 C || Y  (0)  || . 

Hence,  there  exists  V^(t,Y)  such  that 

Ml  <  V L  ( t ,Y)  <  2. 44 Ml  > 

IV^t.Y1)  -  V1(t,Y2)  |  <  2.44||Y1-Y2||, 
Dl4.20)Vl(t,Y)  -  -•5v1(t>Y)- 

In  fact,  for  this  case,  one  such  function  is 


J 

t  1 

1 

1 

j  \ 

II 

/~N 

>< 

' — ' 

> 

4  yt 

"  2 

V  1 

3-/5  | 

l  ] 

!  i 

1  i 

r 

2 

2  | 

/ 

\ 

/ 

and  D^V  (t ,Y)  =  —  V  (Y) .  Solutions  of  (4.21)  satisfv 
1  at  i 

_  ]_ 

iir,''il  =  e  a  || -i '  (0)  !| 

and  we  use 

V2(n')  =  'In'!. 

V2(n')  ,  of  course,  satisfies  a  set  of  conditions  similar  to  (4.23) 

Now  I  can  write 

D(4.19a)VY)  -  --5VY>  +  (2.44)(.5)(/{)V2(n'),  (4.27) 


(4.28) 


(4.19b) 


v2(n’)  <-^v2(n’)  +  ( /|)(/5)(/|)v1(Y) 


Then  using 


I  get 


V3(Y,n  )  =  VX(Y)  + 


u(2.44)(.5)(/|) 


1  -  M  ( •  5) 


v2(n') 


(4.29) 


V1(Y)+ - — u(T5T" V2(n>) 


y(2.44)(.5)(/|) 
1) 


(2.44)  (.5)  (/f)  (/f)  (/5)  (/f) 

+  u - TTTs) - VY) 


<  (-.5  +  -(1^)^-)V3(yV). 


(4.30) 


Hence,  if  uG  [0,  .15],  (4.30)  guarantees  that  (4.19)  is  exponentially  stable. 
Using  linear  time- invariant  system  techniques  one  can  show  that  (4.19)  is 
exponentially  stable  for  u£[0,  1.2).  Thus  we  see  that  the  estimate  of 
range  of  stable  u's  given  by  the  method  of  Theorem  3.9  is  conservative. 

4.3.  Example  Two 

Since  adaptive  control  is  in  general  a  nonlinear  problem  I  consider 
a  nonlinear  example  in  this  section  in  order  to  illustrate  more  fullv  the 
ideas  of  Chapter  3.  I  will,  however,  be  modest  and  choose  the  simplest 
possible  nonlinear  problem.  The  parameters  for  this  simple  example  are 


This  corresponds  to  a  system  for  which  the  high  frequency  gain  of  the  "slow" 

part  of  the  plant  is  known  to  be  y  but  the  desired  value  is  1,  and  hence  I  set 
* 

6^  =  2.  Unlike  the  previous  example  this  plant  is  unstable  without  control 

* 

since  a>0.  The  choice  of  9^  =  -4  provides  the'correct  d.c.  gain  and  the 

k 

steady  state  value  of  e  is  again  zero.  The  choice  of  =  1  is  again  arbitrary 
but,  hopefully,  reasonable. 

This  example  possesses  three  difficulties  which  were  missing  in  the 

* 

first  example.  The  first  difficulty  is  that  because  of  the  feedback  9^,  it  is 
possible  to  make  e  = 0  only  when  xm(0)  =1.  The  second  potential  difficulty 
is  that  the  unperturbed  system  is  nonlinear  and  determining  the  constants 
associated  with  its  exponential  decay  is  not  an  easy  analytic  task.  The 
third  difficulty  is  that  the  region  of  attraction  of  the  exponentially 
stable  zero  solution  of  the  perturbed  system  is  not  the  whole  space  as  it 
was  for  the  linear  example. 

_  * 

I  will  discuss  the  problems  associated  with  the  nonzero  e  first 

since  they  are  the  most  easily  handled  of  the  three  difficulties.  Noting 

k  k 

that  for  this  case  u  =  -4  x  ,  I  can  write  the  tuned  system  as 


(4.32) 


k 

Using  e 


Choosing 


I  get 


°(4. 33)V 


Since  V. 

4 
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x  =  -x  +  1  +  [  .5/  y  0]n 


1 

.* 

r-l 

•5| 

* 

"l" 

yn  = 

1 

L-5 

-1_ 

n  +  m2 

_2_ 

[1  0]n  -  u4/  f 


(x  -  1) 


*  k 

x  -x  ,  I  can  write  the  equations  describing  e  , 


+  [.5/1- 


e  =  -e  +  [  .5/  y  0]n 


(4.33a) 


.* 

un  = 

"-1 

•51 

* 

n  +  y2 

rl 

*  ^  ri 

[1  Q]n  -  y4 /  4  ; 

.  5 

-1 

_2_ 

5  1  2 

(xm  -  1)  .  (4.33b) 


k  & 

x  (0)  =  xm(0)  and  n  (0)  =  0 ,  as  usual,  and  using 


*  *  *  u(.5)(/f) 

V4(e*,n  )  =  | e* |  +  2 


1  -  u(l+  2/5) 


(4.34) 


(e  ,n  )  <  —  ( |  e  |  + 


u(. 5)  (/ b 


1-  u(l+  2/5) 


.  p(.5)(/y)(4)(/|)(/3)  * 

h!>  + - - - - - !e  ' 


1-  n(l+2/5) 


p(.5)(/|)(4)(/|)(/5) 

+ - lx  -  1 1 

1  m  1 


1  -  y(l  +  2/5) 


(4.35) 


.  ,  .  .  d2/5  ,  *  *  u2/5  ,  . , 

<  (-  1+  - — )V  (e  *n  )  + - -  I x  -  1| 

1-  y(l  +  2/5)  ’  1  " 


1  -  y(l  +  2/5) 


*  * 


[e  ,n  )  =  0  at  t = 0  it  follows  from  (4.1),  (4.34),  (4.35)  that 


k  k  k  *-2/5  i  ,  —  f- 

je  I  <  V.  (e  ,n  )<  - -  x  (0)  -  l|e  (4.36) 

-  4  '  1- .,(1+4/5)  m 

for  all  ..  €  [o, - - — ).  Clearly,  if  t  is  sufficiently  small  or  Jx  (0)  -  l| 

1  '  4‘  5  * 

is  sufficiently  ^rall  e  will  be  negligible. 

In  order  to  discuss  the  remaining  aspects  of  this  example  I  need 

T 

the  error  system  for  this  special  case.  Letting  Y  =  *  the  error 

system  can  be  written 


-1 

* 

.  5x 

1  .5Y  Y0 

1  Y  + 

.5/  |  O' 

Y  = 

* 
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■  n* 

(4.37a) 
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In  order  to  simplify  the  problem,  first  consider  the  case  when  x  (0)  =  1. 
k  k 

Then,  e  =  0  and  x  =1.  Under  these  conditions  I  can  bound  u ’  hv 


:  ~  [4. 3||  Y||  +  3.8  ||  Y||  2  + 


Jo 


!IY)I 


-V  J 


JO 


(4.39) 


Now,  I  consider  the  fast  and  slow  subsystems  separately  in  order  to  find  a 
Lyapunov  function  for  each.  Since  the  fast  subsystem  is  again  given  by 


(A. 21),  I  will  use  l^Cn')  for  this  case  also.  The  slow  subsystem  is  given 
by 


(A.  4 


The  zero  solution  of  this  nonlinear  slow  subsystem  is  exponentially  stable 
in  the  large.  However,  it  is  not  easy  to  analytically  determine  the 
constants  associated  with  the  exponential  convergence  of  (A. AO).  Therefore, 
I  resorted  to  numerical  techniques.  Since  I  wanted  =  (1,0)  ° 

(x'.ep  =  (0,A)  to  be  in  the  region  of  attraction  of  the  zero  solution  of 
(A. 37),  I  choose  to  estimate  the  exponential  rate  of  convergence  of  (A. AO) 
for  ||Yo||  <  5.  First,  I  simulated  (A. AO)  with  initial  conditions  at  16 
uniformly  spaced  points  on  a  circle  of  radius  5  centered  at  Y  =  0.  Second, 

I  plotted  the  maximum  value  of  the  |Y|  over  the  16  trajectories  at  each 
sample  instant.  The  integrations  were  done  using  the  IMSL  subroutine  DVERK 
which  determines  its  own  step  size  to  remain  within  a  prescribed  tolerance 
in  integration  error.  For  plotting  purposes  I  sampled  the  integration  10 
times  per  second.  Finally,  I  used  the  plot  of  the  maximum  value  of  jj Y| 
versus  time  to  estimate  that 


jj Yj|  <  1.82  e  °-077t  J Y ( 0 ) J  ,  VY(0)GS5  ,  (A. 

where  S<-  =  {Y:  j|Y||  <  5}.  Figure  A.  3  shows  the  plot  of  the  maximum  j|  Y||  and 

the  bound  (A.A1).  Figure  A. A  shows  j| Yj|  vs  t  and  Y^  vs  Y^  for  an  example 

T 

trajectory  with  initial  condition  Y(0)  =  [0,5]  . 
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Figure  4.4.  (a)  Plot  of  vs  Y7  and  (b)  plot  of 

trajectory  of  (4.40)  beginning  at  Y=  [0,5] 


vs 

alT 


From  (4.41)  and  Theorem  3.7,  it  follows  that  there  exists 
V^(t,Y)  such  that 

1Y!|  <  V5(t,Y)  <  1.82  ||Y|j  ,  VYGS5  , 

|V5(t,Y1)  -  V5(t,Y2)|  <  L(t  ,5)|  Y1  -  Y2||  ,  Y16S5,Y26S5  ,  (4.42) 

D_(4.40)V5(t’Y)  1  -q(.077)V5(t,Y)  ,  0<q<l  . 

For  the  purposes  of  this  discussion  it  will  suffice  to  make  the  unjustified 
(and  optimistic)  assumptions  that 

q  =  1  ,  L(t ,5)  <  1.82  .  (4.43) 

Using  (4.39),  (4.42),  and  (4.43)  I  write 

°"(3.37a)e*  =  oV5(t,Y)  1  "  -077  Vt,Y)  +  1-82  ( -5)  (/|)V2(n ' )  ,  (4.44) 

°(3.37b)eSoV2(,1')  i(-±+2/S+  (,3)/5  |Y|)»2(-,’) 

(4.45) 

+  (4.3  +  3.8  |  Yjj  +  v  2  | Y||  2) V-  (t ,Y)  . 

Letting  be  an  upper  bound  on  ||Y||,  I  define 

(.91)/  4  (4.3  +  3.8  M  +  y'i  M2) 

=  .077  -  u  - 1 - - - 2 - £  .  (4.46) 

°  1  -  u(2/5  +  .077  +  .5/5  M  ) 

o 

It  is  easily  shown  that  if  Mq  is  an  upper  bound  on  J Yj|  ,  then 

(•91)/| 

V  (t,Y,n')  =  V_(t,Y)  +  u  - - - - - — -  V?  (n  ' )  (4.47) 

3  1  -  y  (2/5  +  .077  +  .5/5  M  )  ~ 


satisfies 


°(4.37)e*  =  oV6^t,Y,n'^  -  a6V6  ^  ,Y  ’ n  ' ) 


(4.48 


Hence,  if  a  >0  and  M  is  an  upper  bound  on  V, (0 ,Y(0) ,n ' (0) ) ,  then  M 

DO  D  O 

is  an  upper  bound  on  ||Y||  since  j|  Y||  <  V,.(t,Y)  <  V^t.Y.q')  -  V&  (0 ,  Y  (0)  ,  n  '  (0) ) 

This  is  the  motivating  force  behind  the  definition  of  M(Y , • ')  as  the  solution 


M(Y  ,n  ')  =  1.82  H  Y||  + 


u(.91)/|  |n'| 

1  -  u(2/5  +  .077  +  .5/5  M(Y.n') 


(4.49 


Then  taking  Mq  =  M(Y(0) , n ' (0) )  ensures  that  Mq  is  an  upper  bound  on 
(0 , Y(0) , n ' (0) ) .  If  I  want  to  guarantee  by  this  method  that  the  set 
(Y,  n ' )  €  {Y,  n  '  :  Y€s„,nf  =  0}  is  in  the  region  of  attraction  of  the  zero 
solution  of  (4.37),  I  must  show  >  0  for  =  9.1.  This  can  be  done  only 
for  u€  [0,  .00034] . 

Returning  to  the  more  general  problem  when  xm(0)  ^  1  but  for 

example  satisfies  |xm(0)  -  l|  <2,  one  expects  that  the  requirements  on 

u  will  be  at  least  as  restrictive  as  in  the  x  (0)  =  1  case.  Thus,  start  with 

ra 

1  *  1 

u  <  .00034.  Clearly  |e  j  is  then  negligible  with  respect  to  the  other  terms 
in  (4.38)  so  that  the  bound  (4.39)  is  good  with  respect  to  e  at  least. 

t  *  I 

Because  | e  |  is  negligible 


x  =  x  =  1  +  (x  (0)  -  l)e  C 
m  m 


(4.50 


In  order  to  truly  determine  the  effects  of  this  time  varying  x  one  may 
need  to  find  a  new  bound  similar  to  (4.41)  for  the  slow  subsystem 
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A 

—  ” 

-1  .5x 

•5Y1Y2 

* 

Y  + 

2 

_-x  0 

-  -Yl- 

(4.51) 


with  x  given  by  (4.50).  However,  I  know  that  the  zero  solution  of  (4.51) 

*  * 

is  exponentially  stable  in  the  large,  provided  x  ?0.  Furthermore,  x  -  1  for 

t >  3.  Thus,  for  t>  3  the  previous  analysis  should  hold  with  Mq  =  M(Y(3), 

n ' (3) ) .  With  these  thoughts  in  mind  one  might  consider  a  pseudo  worst  case 

analysis  in  which  one  assumes  that  V^(t,Y)  satisfies  only  V5  (t  ,Y)  <  0  for 

0  <  t  <  3.  Under  this  assumption  and  the  assumption  that  (4.39)  is  a 

reasonable  bound  over  0  <  t  <  3, one  gets  j.,  >  -  .077  for  p£[0, .00034]  and 

—  -  0 

M  =9.1.  From  this  it  follows  that 
o 

V(3 ,Y ,n ' )  <  e(,077) (3)V(0,Y,n’)  . 


Since  M  =  9.1  must  bound  V(3,Y,n'),  this  psuedo  worst  case  analysis  reveals 

that  the  guaranteed  region  of  attraction  is  reduced  by  a  factor  of 
-.231 

e  =  0.8. 

How  useful  is  it  to  know  that  for  u  <  .00034  I  can  guarantee  that 
my  desired  region  is  contained  in  the  region  of  attraction?  This  requires 
the  unmodelled  high  frequency  dynamics  to  be  3,000  times  as  fast  as  the 
modelled  dynamics.  This  seems  like  a  pretty  unrealistic  requirement. 
However,  to  be  accurate,  the  answer  to  the  first  question  depends  on  the 
answer  to:  "How  big  can  u  really  be  before  the  region  of  attraction  does 
not  contain  my  desired  region?"  Working  with  the  desired  region 
(Y,i) €  (Y,n  :YG  S_ ,n  =  0} ,  I  tested  via  simulation  whether  u  =  .1  was  small 


enough  so  that  the  desired  region  was  contained  in  the  region  of  attraction 


of  the  zero  solution  for  the  x  (0)  -  1  case.  I  found  that  it  was  small 

m 

*  & 

enough.  For  this  test  I  integrated  (4.37)  with  x  =  1  and  e  =  0  using  the 

same  16  initial  conditions  on  Y  as  when  I  studied  (4.40)  and  the  initial 

condition  r,'(0)  =  0.  Then,  I  again  took  the  maximum  of  ||Y||  at  each  sample 

instant.  Figure  4.5  shows  the  plot  of  max  |Y|  vs  t  for  both  the  unperturbed 

system  (4.40)  and  the  perturbed  system  (4.37).  Figure  4.6  shows  a  sample 

T 

trajectory  of  the  perturbed  system  with  initial  conditions  Y (0)  =  [0,5]  and 
rf(0)  =  0  and  with  p  =  0.1.  Compaing  0.1  and  0.00034  it  is  obvious  that  the 
estimate  provided  by  the  proof  of  Theorem  3.9  is  so  conservative  that  it  is 
of  little  or  no  practical  use. 

4.4.  Discussion 

I  have  investigated  two  simple  examples  which  have  illustrated  the 
theory  of  Chapter  3.  That  is,  for  each  example  there  was  a  p  >  0  such  that 

-  it 

for  iiS  [0,u)  the  zero  solution  of  the  error  system  with  e  =  0  retained  its 
exponential  stability.  Furthermore,  it  was  shown  that  p  does  not  have  to 
be  an  incredibly  small  number.  However,  the  examples  also  showed  that,  even 
for  these  simple  systems,  the  estimates  of  p  and  the  region  of  attraction 
developed  from  the  techniques  of  the  proofs  of  the  theorems  in  Chapter  3  are 
so  conservative  that  they  are  of  little  practical  use.  This  conservation 
arises  from  the  fact  that  the  technique  assumes  the  terms  coupling  the  fast 
and  slow  subsystems  together  always  have  the  maximum  detrimental  effect  upon 
the  stability  of  the  error  system  while  the  coupling  terms  actually  have  a 
much  smaller  effect  on  the  stability  of  the  error  system. 


As  the  system  gets  larger  and  contains  more  adjustable  parameters, 
this  conservativeness  problem  will  only  get  worse.  This  is  because  of  the 
time-varying  nature  of  the  larger  systems.  First,  the  technique  bounds  a 
time-varying  quantity  with  its  maximum  absolute  value.  Then,  it  applies 
the  maximum  detrimental  effect  idea.  This  double  maximum  approach  can 
hurt  the  accuracy  of  the  estimate  very  much  since  the  effect  of  time-varying 
coupling  term  may  depend  upon  its  average  value  or  the  average  value  of  its 
correlation  with  another  time-varying  quantity. 


CHAPTER  5 


CONCLUSION 

This  study  of  model  reference  adaptive  control  in  the  presence 
of  unmodelled  high  frequency  dynamics  has  shown  several  interesting  results. 
The  first  of  these  results  is  that  for  any  input  in  the  class  of  signals 
which  provide  exact  sufficient  richness,  the  system  has  a  unique  e  =  0 
equilibrium  provided  the  high  frequency  dynamics  are  at  frequencies  above 
the  frequencies  in  the  input.  Along  with  this  result,  it  was  shown  that 
the  equilibrium  values  of  the  parameters  are  only  a  small  distance  from 
the  ideal  parameters  for  the  unperturbed  case. 

The  remaining  results  concern  the  stability  of  this  equilibrium 

Taking  the  two-time  scale  approach  to  the  problem  leads  to  the  satisfying 

theoretical  result  that  this  equilibrium  is  exponentially  stable  after  the 

model  transients  have  settled  provided  that  the  high  frequency  unmodelled 

dynamics  are  sufficiently  high  frequency.  Through  the  remarks  it  was 

pointed  out  that  this  technique  can  be  used  to  show  that  the  error  system 

* 

converges  exponentially  to  an  0(a)  region  around  e  =  0,  9  =  9  provided 
u  is  sufficiently  small. 

Finally,  the  examples  showed  that  there  indeed  exists  a 
a  >0  such  that  the  equilibrium  is  exponentially  stable  for  all  y 6  [0,y). 
However,  the  examples  also  showed  that  the  estimates  available  from  the 
proofs  are  not  large  enough  to  be  of  practical  use.  Clearly,  an  important 
goal  for  new  research  is  to  find  an  approach  to  this  problem  which  can 
result  in  better  estimates  of  the  range  of  u  for  which  the  equilibrium 
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retains  its  stability.  Hopefully,  an  approach  which  provides  realistic 
estimates  will  come  close  to  providing  a  fundamental  understanding  of 
the  problem  and  thus  lead  to  measures  for  counteracting  the  destabilizing 
effects  of  the  unmodelled  high  frequency  dynamics. 


& 


■J 


APPENDIX 


PROOFS  OF  THEOREMS  3.6  AND  3.7 


The  following  proofs  are  taken  directly  from  Yoshizawa  (1966) 
Proof  of  Theorem  3.6.  Let  V(t,x)  be  defined  by 

V(t,x)  =  sup  I  x(t  +  t  ;x,t)jj  eCT  .  (A. 

x>0 

Then  clearly  |j  x|  <  V(t,x)  and  by  (3.7) 

V(t,x)  <  sup  Ke  CT|x|eCT  =  K||  xj[  .  (A. 

x>0 

Since  the  system  is  linear,  we  have  the  relation 
x(t+x;x,t)  -  x(t  +  x;x',t)  =  x(t  +  x  ;x  -  x' ,t)  ,  and  hence, 

j  V(t  ,x)  -  V(  t  ,x '  )  |  <  sup  |x(t+t;x,t)  -  x(t  +  x;x' ,t)j|eCT 

r>0 

<  sup  Ke  CT||  x-x' |  eCT  =  K||  x  -  x '  |  .  (A. 

T>0  '  ’’ 

Now  we  shall  prove  the  continuity  of  V(t,x).  Take  a  6>0. 

We  have 


|V(t  +  '5,x')  -  V(t,x)l  <  |v(t  +  5,x')  -  V(t+5,x)| 

+  jv(t  +  <5,x)  -  V(t  +  6  ,x(  t  +  6  ;x,  t) )  | 

+  ;V(t + 5 ,x(t + 6 ;x,t))  -  V(t ,x) |  .  (A. 


Since  V(t,x)  is  Lipschitzian  in  x  and  x(t  +  <5;x,t)  is  continuous  in  $ ,  the 
first  two  terms  are  small  when  II  x  —  x '  f|  and  5  are  small. 


Let  us  consider  the  third  term.  Since  x(t + c + x ;x(t +  5 ;x,t) ,t + 6) 
=  x(t  +  <5  +  t;x,t),  we  have 


|  V(t  +  5  ,x(t  +  <5  ;x,  t)  )  -  V(t,x)|  =  jsup  ||  x(t  +  6  +  x  ;x(t  +  5  ;x,  t)  ,  t  +  6  )  ||  eC 

T>° 

-  sup  |  x ( t  +  t  ; x , t ) ||  eC‘  j 

T>0 

=  jsup  |  x(t  +  T  ;x,t)||  eCTe  C' 

t>Q 

-  sup  |x(t+  r;x,t)J|eCT|  .  (A. 5) 

T>0 

Set  a(5)  =  sup  ||  x(t  +  x  ;x,  t) ||  eCT  .  Then  a(5)  is  nonincreasing  and  a(£)-*0 
T>0 

as  S-’-O,  because  |  x(t  +  T  ;x,t)||  e  '  is  a  bounded  continuous  function  for  all 
T  >  0 .  Thus  , 

j  V(t  +  <5  ,x(t  +  5  ;x,  t) )  -  V(t,x)|='a(i)eCJ-a(0)|  (A. 6) 


implies  that  the  third  terra  on  the  right-hand  side  of  (A. 4)  tends  to  zero 
as  o-*0.  Therefore,  the  continuity  of  V(t,x)  is  verified. 

Finally,  we  shall  establish  condition  (iii) .  Let  x'  =  x(t +  h;x,t) 

h  > 0 .  Then, 


V  ( t  +  h  ,  x ' )  =  sup  ||x(t  +  h  +  t  ;x' ,t  +  h)j|  eC  ‘ 

t>0 

=  sup  ||x(t  +  r;x,t)|ec  e  Ch  <  V(t,x)e  Ch, 

r>h 

which  implies 

vfeHa'UV);.,)  ,  v(  ,  Spl 

h  -  h 

From  this,  we  obtain  D*  _.(t,x)  <  -  cV(t,x). 


(A.  7) 


(A. 8) 


Proof  of  Theorem  3.7.  For  a  constant  q  such  that  0<q<  1,  define 


V(t,x)  =  sup  |x(t  +  -r;x,t)||eqaT 
T  >  0 


(A. 9) 


I  a 

Then  |J  x|j  <  V(t,x)  and  for  x£Sa  =  {x:||x|  <  6} 


xj  <  V(t,x)  <  sup  K(B)e~(1  q)aT|xj|  <  K(B)|: 

T>0 


(A. 10) 


Now  we  shall  prove  (ii) .  In  the  course  of  the  proof,  we  shall 
determine  L(t,S)  explicitly.  Let  T(6)  be  such  that  K(g)  =  e^  9)aT(£0. 

If  t>T(8),  then  K(S)e  ^  q^ai||x||  <  [|  xj|  and  hence,  from  (A. 10),  it  follows 

that  V(t,x)  must  be  defined  for  i  such  that  0  <  t  <  T(ct)  .  Therefore,  for 

x£S  and  x'6S„, 

P  p 

:V(t,x)  -  V(t,x')(<  sup  |jx(t  +  x;x,t)  -  x(t  +  T;x',t)||  eqaT  ' 

0  <  t  <  T(S)  " 


< 


gqaT(3) 


exp  [ 


t  +  T(6) 

/  M(s,B)ds]j|x 

t 


(A. 11) 


where  M(s , 8) 

•  M  ( s  ,  3  )  |  x  - 


is  a  continuous  function  of  s  such  that  jjf(s,x)  -  f(s,x') 

*’!  for  .SS(K(!)  and  >'£S1K([). 

The  remainder  of  the  proof  can  be  verified  by  the  same  * 


argument  as  in  the  preceding  proof. 
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